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EXECUTIVE  SUMMARY 


STRESS  WAVE  PROPAGATION  IN  UNSATURATED  SANDS:  CENTRIFUGE  MODELING 

Explosive  model  testing  wae  conducted  using  a  centrifuge  in  order  to 
simulate  prototype  stresses  and  ground  motions  in  a  representative 
cohesionless  backfill.  Models  viere  constructed  of  sand  and  compacted 
moist  to  a  constant  void  ratio  using  a  vibratory  technique.  Exploding 
detonators  were  used  to  simulate  contained  bombs  in  the  backfill  material. 
The  objective  of  this  study  wae  to  determine  the  influence  of  moisture 
content,  at  the  time  of  backfill  compaction,  on  blast  induced  stress  wave 
propagation.  Models  were  constructed  to  1/18.9  and  1/26.3  scales. 
Explosives  consisted  of  1031  mg  and  350  mg  of  PBX  9407  and  were  buried  to 
depths  of  7.6  cm  and  5.4  cm  respectively.  Accelerating  these  scaled 
models  to  18.9  g's  and  26.3  g's  simulated  prototype  charges  of  6.9  kg  and 
6.4  kg  (7.8  kg  and  7.3  kg  TNT  equivalent)  at  a  depth  of  burial  of  1.4 
meters.  peak  stress,  peak  particle  acceleration,  and  peak  particle 
velocity  are  presented  as  a  function  of  scaled  distance  for  both  Tyndall 
beach  and  Ottawa  20-30  sands.  Attenuation  coefficients  (n)  for  this  data 
show  some  influence  with  respect  to  saturation.  Peak  stress  and  peak 
particle  acceleration  intercepts  (b)  at  a  scaled  distance  of  one  (R/W'^  > 
1  m/kg''’  and  1  ft/lb''’)  for  Tyndall  sand  are  lowest  at  0  and  53  percent 
saturations  and  are  maximum  at  35  percent  saturation.  Peak  particle 
velocity  intercepts  are  lowest  at  0  and  70  percent  saturations  and  again 
are  maximum  at  35  percent  saturation.  Peak  stress  Intercepts  for  Ottawa 
sand  follow  a  similar  trend  to  stress  intercepts  of  Tyndall  sand. 
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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

The  general  objective  of  tMa  reaearch  program  is  to  determine  the 
influence  that  the  degree  of  saturation  during  compaction  of  aand  has  on 
blast-induced  ground  shock  and  stress  wave  propagation.  Pour  specific 
objectives  arise  out  of  this  general  objective. 

1.  Develop  empirical  equations  that  accurately  predict  peak 
stresses,  peak  particle  velocities,  and  peak  particle 
accelerations  in  Tyndall  beach  and  Ottawa  20-30  sands  as  a 
function  of  compactive  saturation. 

2.  Compare  centrifuge  results  with  work  on  dynamically  loaded 
sands  from  the  split-Hopkinson  pressure  bar  reported  by  Ross 
et  al.  <1988),  Pierce  et  al.  (1989)  and  Charlie  et  al.  (1990). 

3.  Compare  the  validity  of  centrifuge  results  with  those  of  full 
scale  blasting  by  Drake  and  Little  (1983). 

4.  Compare  centrifuge  results  with  unsaturated  soil  mechanic 
theory  proposed  by  Fredlund  (1985). 

B .  BACKGROUND 

Blast  loading  of  soils  is  of  considerable  importance  to  those 
concerned  with  buried  structures  because  of  the  high  intensity,  high 
amplitude  compressive  stress  waves  produced.  Compacted  backfill,  placed 
around  and  over  a  buried  structure,  transmits  the  blast  induced  stresnes 
and  thus  determines  the  survivability  and  vulnerability  of  the  structure. 

Althou^jh  Drake  and  Little  (19B3)  conducted  numerous  blast  tests  and 
developed  stress  and  ground  notion  prediction  equations  on  various  soils, 
there  has  not  bean  a  systematic  explosive  testing  program  conducted  on 
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unsaturated  sands.  Laboratory  research  on  sands  by  Rosa  st  al.  (1988), 
Pierce  et  al.  (1989),  and  Charlie  et  al.  (1990),  utilising  the  split- 
Kopkinson  pressure  bar,  has  shown  that  the  degree  of  saturation  during 
compaction  is  influential  on  dyr.amic  stress  wave  velocity  and  stress 
transmission  ratio. 

C.  SCOPE 


This  report  presents  the  results  of  systematic  explosive  testing  on 
unsaturated  Tyndall  beach  and  Ottawa  20*30  sands.  Our  purpose  is  to 
determine  the  influence  of  moisture  content  during  compaction  on  blast 
induced  stress  and  ground  motion  response.  Explosive  events  are  carried 
out  in  small-scale  models.  A  centrifuge  is  utilized  to  subject  the  models 
to  increased  accelerations  such  that  prototype  blast  induced  stresses  can 
be  simulated.  Scaling  laws  are  used  to  relate  model  and  prototype 
performance. 

performing  scaled  model  testing  on  the  centrifuge  allows  prototype 
explosive  events  to  be  conducted  at  a  fraction  of  the  time  and  money 
required  to  perform  large  scale  blast  testing.  Charge  masses  of  1031  mg 
and  350  mg  of  PBX  9407  are  utilized  to  induce  explosive  loadings  in  models 
of  1/18.9  and  1/26.3  scale.  Accelerating  the  models  to  18.9  and  26.3 
times  the  acceleration  of  earth's  gravity  results  in  scaled  prototype 
charges  of  6.9  kg  and  6.4  kg  respectively  (7.B  kg  and  7.3  kg  TNT 
equivalent)  detonated  at  a  depth  of  burial  of  1.4  meters. 
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SECTION  II 


LITERATURE  REVIEW 

This  chapter  briafly  provides  insight  on  the  theories  and  principles 
involved  in  this  research  effort.  These  theories  and  principles  are  the 
basis  for  drawing  conclusions  in  this  report. 

A.  STRESS  WAVES  AND  GROUND  MOTIONS 


1.  Stress  Wave  Propagation 

As  a  stress  wave  travels  further  from  its  source,  the  wave 
front  becomes  less  curved  and  particles  traveling  in  a  direction  away  from 
the  source  are  assumed  to  os  parallel.  Representing  wave  motions  in  this 
manor  is  known  as  one-dimansional,  or  planar  wavs  propagation.  Planar 
wave  propagation  and  constrained  modulus,  M,  can  be  related  by  considering 
a  stress  wave  propagating  down  a  bar.  Derivation  begins  by  using  Newtons' 
Second  Law  of  motion 

F, -jna,  (2.1) 

where  F.  represents  an  induced  force  on  an  element  of  the  bar  of  mass  m, 
and  a,  is  the  acceleration  of  the  element  (Dowding,  1965).  Applying  this 
force  over  an  elemental  area  of  the  bar  results  in  the  equation 

o  •  p  c,,  *  (2.2) 

where  a  is  stress,  p  is  the  density  of  the  material  in  the  bar,  c,  is  the 
plana  wavs  propagation  velocity  through  the  bar  (or  seismic  velocity),  and  it 

is  the  velocity  of  the  element,  or  particle  velocity  (Kolsky,  1963).  The 
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quantity  pc,  in  Equation  (2.2)  is  coitanonly  known  as  the  acoustic  impedance 
of  the  medium.  Equation  (2.2)  predicts  stress  on  the  stress  wave  front  of 
planar,  spherical  and  cylindrical  stress  waves  (Rinehart,  197S). 

Considering  elastic  materials,  strain  can  be  expressed  in  terms  of 
particle  velocity  and  seismic  velocity 


c  = 


(2.3) 


Seismic  velocity  and  constrained  modulus  can  be  related  in  the  form 


Similarly,  shear  wave  velocity  is  expressed  in  terms  of  shear  modulus,  G„ 


(Rinehart,  1975).  Shear  modulus  and  constrained  modulus  are  related  to 
elastic  modulus  and  Poisson's  ratio,  with 


E 


2(1  ♦  p) 


and 


(2.6) 


M 


g(l  -  U) 

(1  ♦  p)  (1  -  2\^) 


(2.7) 


(Lambe  end  Whitman,  1969). 

Seismic  velocity  is  measured  at  strains  below  10*^  percent  and  may  be 
used  as  a  crude  index  for  predicting  stress  wave  propagation.  Generally, 
low  seismic  velocity  is  associated  with  low  relative  density  and  indicates 
poor  stress  wave  transmission.  Have  attenuation  can  be  described  as  the 
change  in  peak  stress,  peak  particle  velocity,  or  peak  particle 
acceleration  with  time  or  distance  and  is  greatest  among  low  density 
soils.  Attenuation  coefficients,  n,  are  determined  from  the  slopes  of 
least  squared  regression  lines  of  blast  data  when  plotted  on  log-log 
scale.  Table  2.1  shows  gror.nd  shock  parameters  for  sands  of  various 
saturation  and  relative  density  and  determined  from  several  hundred 
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exploaive  tests  by  Drake  and  Little  (1983).  Seismic  velocities  are  shown 
in  column  five  (5)  of  Table  2.1.  and  acoustic  impedance  and  attenuation 
coefficients  are  shown  in  columns  six  (6)  and  seven  (7). 

2.  Charge  Mass  and  Distance  Relationship 

To  analyze  blasting  data  when  charge  mass,  W,  and  distance,  R, 
from  the  blast  vary,  a  scaling  technique  relating  the  two  variables  is 
useful.  Two  common  approaches  are  the  cubed  root  scaling,  R/w'^,  and  the 
square  root  scaling,  (Dowding,  1985).  Cubed  root  scaling  is  based 
on  explosive  charge  distribution  of  spherical  form,  and  square  root 
scaling  assumes  a  cylindrical  explosive.  Throughout  this  report  the  cubed 
root  scaling  technique  will  be  used. 

3.  Ground  Shock  Coupling  Factor 

The  magnitude  of  blast-induced  stresses  and  ground  motions  is 
greatly  enhanced  as  a  weapon  penetrates  more  deeply  into  the  soil  before 
it  detonates.  The  ground  shock  coupling  factor,  f,  is  a  relationship  of 
blast  energy  transmitted  to  the  surrounding  medium.  It  is  defined  as  the 
ratio  of  the  ground  shock  magnitude  from  a  partially  buried  or  shallow 
weapon,  to  that  from  a  fully  buried  weapon  in  the  same  medium.  Figure  2.1 
shows  ground  shock  coupling  factors  as  a  function  of  scaled  depth  for 
blasts  in  air,  soil,  and  concrete  (Drake  and  Little,  1983).  Depth  of 
burial  is  meas\.  i  center  of  mass  of  the  expulsive. 

4.  Ground  Shock  Prediction  Equations 

Free-field  stresses  and  ground  motions  can  be  characterized  by 
rapidly  decaying  exponential  time  histories  that  decay  monotonically,  or 
attenuate,  as  they  propagate  outward  from  the  explosion.  These 
attenuating  ground  motions  can  be  represented  by  a  straight  line  when  the 
data  is  plotted  on  log-log  scale.  The  equation  for  a  line  on  log-log 
scale  has  the  form 


6 


OMOUND  IHOCk  COUI>LINO  FACTOR 


logy  ■  SI  logx  logh 


(2.8) 


where  y  !■  the  ordinate  (i.e.  peak  streae,  peak  particle  acceleration,  or 
peak  particle  velocity),  and  x  ia  the  abacisaa  or  scaled  distance  from  the 
charge.  The  constants  m  and  b  represent  the  slope  of  the  line  and  y- 
intercept  of  the  line  (at  a  scaled  distance  of  one  on  log-log  scale) 
respectively.  The  inverse  log  of  Equation  2.8  results  in  an  equation  of 
the  form 

y  ■  b  x".  (2.9) 

Pree-field  etreae  and  ground  motion  data  presented  by  Drake  and 
Little  (1983)  is  in  the  form  of  Equation  (2.9).  Equations  for  peak  stress 
(P),  peak  particle  velocity  (V),  peak  particle  acceleration  (a),  peak 
particle  displacement  (d),  and  impulse  (I),  are  as  follows: 


P  »  f  0.049  pCe|2. 52-^1  (kPa) 

J 


y  =  f  160 


f-^r 


(ft/8) 


(2.10a) 


(2.10b) 


(2.11a) 


V  =  f  49 


(m/B) 


(2.11b) 


K’/*  =  f  50  c. 


(g'a-ENGLlSH) 


a  K’-'*  =  f  126  12.52 


(g's-si) 


(2.12a) 


(2.12b) 


8 


(ft) 


(2.13a) 


d 

ifi/i 


=  f  500 


1 

c 


d 


=  f  60  — 
c 


(m) 


I 

^i/j 


=  f  p  1.1 


(psi-s) 


=  fp0.19(2.52^)'’'  (kPa-s) 

where 

a  ■  peak  particle  acceleration  (g's); 

Cj  *  seiemic  velocity  (ft/e;  m/e); 
d  •  peak  particle  displacement  (ft;  m) ; 
f  «  ground  shock  coupling  factor; 

g  •  earth  gravity  force  (1  g  *  32.2  ft/s’  "  9.81  ro/e^); 

I  •>  peak  impulse  (psi-s;  kPa-s); 
n  «  attenuation  coefficient; 

P  «  peak  stress  (psi;  kPa); 

R  *  distance  to  the  explosive  (ft;  m) ; 

p  ■  mass  density  (Ib-eVft^;  kg/m’); 

pCj  •  acoustic  impedance  (Ib-e/ft’;  kg/m^-s); 

V  ■  peak  particle  velocity  (ft/s;  m/s);  and 
W  *  charge  weight;  charge  mass  (lb;  kg). 

e 

By  setting  pc.  to  units  of  Ib-s/in’-ft,  Equations  (2.10a)  and  ( 
be  combined  with  the  common  term  f  160 (R/H''’)'"  to  yield 

P  *  p  V. 

Equation  2.15  is  the  same  as  Equation  (2.2)  shown  earlier. 


(2.13b) 


(2.14a) 


(2.14b) 


2.11a)  can 


(2.15) 
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B. 


UNSATURATED  SANDS 


1.  StresB  Wave  Propagation  in  Unsaturated  Sande 

Currently,  there  are  no  theoretical,  empirical  or  numerical 
methods  that  predict  stress  or  ground  motion  behavior  in  cohesionlesa  soil 
as  a  function  of  saturation.  Most  field  specifications  for  cohesionless 
soil  placement  is  based  on  the  assumption  that  saturation  has  little  or  no 
influence  on  the  dynamic  behavior  of  the  soil. 

Robb  et  al.  (1988)  and  Pierce  et  al.  (1989)  performed  high-amplitude 
split-Hopkinson  pressure  bar  (SHP6)  tests  on  Eglin  and  Ottawa  20-30  sands. 
They  determined  that  moisture  content  and  capillary  pressure  during 
compaction  and  testing  were  important  in  predicting  wave  velocity  and 
stress  wave  transmission.  The  results  of  Ross  et  al.  (1988)  showed  that, 
for  sands  compacted  moist  to  a  constant  dry  density,  wave  velocity  and 
stress  transmission  increased  as  the  saturation  was  increased  from  0  to  20 
percent,  remained  steady  to  70  percent  saturation,  and  then  decreased 
(Figures  2.2(a),  (b),  (c),  and  (d) } .  Pierce  et  al.  (1989)  conducted 
similar  tests  to  those  of  Ross  et  al.  (1988)  except  that  the  specimens 
were  compacted  both  moist  to  a  constant  dry  density,  and  compacted  dry  to 
a  constant  density,  then  wetted.  Pierce  et  al.  (1989)  results  for  Eglin 
sand  compacted  moist  (Figures  2.2(b)  and  (d))  are  similar  to  those 
obtained  by  Ross  et  al.  (1988),  and  results  for  specimens  compacted  dry 
(Figures  2.2(e),  (f),  (g),  and  (h))  indicate  little  change  in  wave 
velocity  and  stress  transmission  with  saturation. 

Similar  results  ware  reported  by  Charlie  at  al.  (1990)  from  SHPB 
teste  conducted  on  silica  50/80  sand  compacted  moist  to  a  constant  dry 
density.  Transmitted  stresses  increased  from  0  to  approximately  20 
percent  saturation,  then  decreased  with  increasing  saturation  (Figure 
2.3(a)).  Have  velocities,  however,  remained  steady  to  approximately  70 
percent  saturation,  than  decreased  with  increased  saturation  (Figure 
2.3(b)).  They  concluded  that  both  water  content  and  dry  density  need  to 
be  controlled  during  compaction  if  stress  transmission  and  attenuation  are 
critical. 
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Figure  2.2  Resulta  obtained  by  Roee  et  al.  (1988)  and  by  Pierce  et  al.  (1989)  on  Eglin  and  Ottawa  aanda 
compacted  molat.  (a)  Ottawa  wave  speed.  (b)  Eglin  wave  epeed.  (c)  Ottawa  tranemiaslon 
ratio.  (d)  Eglin  trananiisaion  ratio.  Reeulta  of  Pierce  et  al.  (1989)  taken  front  Charlie 
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Charlie  and  Walah  (1990)  conducted  one-dimensional  (planar)  blast 
testing  in  soil  models  on  the  centrifuge  simulating  full  scale  blast 
events.  Wave  velocity  results  showed  an  increase  from  0  to  20  percent 
saturation,  remained  steady  to  60  percent  saturation,  and  then  decreased 
(Figure  2.4).  Unpublished  results  of  average  pea)t  stress  for  various 
scaled  depths  and  explosive  mass  are  shovm  in  Figure  2.5.  The  general 
trend  is  an  increase  in  peak  stress  as  saturation  increases  from  0  to  40 
or  60  percent,  then  peak  stresses  decrease  as  saturation  increases  to  80 
percent. 

Veyera  and  Fitzpatrick  (1990)  examined  the  relationship  between 
compaction  moisture  content,  soil  microstructure  and  dynamic  stress 
transmission  from  tests  conducted  on  the  SHPB  fcr  Ottawa  20-30  sand.  They 
found  that  the  compactive  energy  required  to  obtain  constant  dry  density 
is  strongly  dependent  on  the  amount  of  moisture  present  during  compaction 
(Figure  2.6).  They  attribute  this  to  the  formation  of  preferred  particle 
orientation  and  capillary  pressure  during  compaction. 

Dynamic  shear  modulus,  G.,  is  a  principal  soil  property  required  for 
evaluation  of  wave  propagation  in  fins  grained  cohesionless  soils.  Wu  at 
al.  (1964)  determined  the  significance  of  capillarity  effects  on  shear 
modulus  and  the  effective  grain  size  D,o.  Results  showed  a  significant 
increase  in  shear  modulus  when  the  soils  were  compacted  moist  to  a 
constant  dry  density.  These  effects  were  greatest  among  soils  having  the 
smallest  grain  diameter.  Resonant  column  teste  or.  glacier  way  silt 
(Figure  2.7)  show  that  shear  modulus  displays  a  definite  peak  between  10 
and  20  percent  saturation  and  for  various  confining  pressures.  This 
increase  in  G,  is  a  result  of  an  increase  in  shear  wave  velocity. 

2.  Effective  Stresses  in  Unsaturated  Sands 

Stress  transmission  is  inherently  related  to  interparticle 
contact  and  normal  contact  stress.  Shukla  and  Prakash  (1988)  have 
conducted  blast-induced  stress  wave  propagation  tests  in  plate  specimens 
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and  have  shown  that  wave  speed  drops  with  increasing  porosity.  They  go  on 
to  conclude  that  wave  velocity  shows  a  strong  dependence  on  the 
microstructure  of  the  porous  media. 

Stress  waves  are  transmitted  through  the  soil  skeleton  across  the 
small  mineral  contacts.  The  response  of  a  eoil  mass  to  changes  in 
compressive  stress  depends  a  great  deal  on  the  intergranular  stresses. 

The  equation  for  effective  stress,  developed  by  Terzaghi  (1943),  has 
the  form 

o'  ^  a  ~  u  (2.16) 

where  o'  denotes  effective  stress,  a  is  the  normal  stress,  and  u  is  the 
pore  fluid  pressure.  In  the  case  of  partially  saturated  soil  where  pore 
spaces  contain  two  fluids,  (water  and  air).  Bishop  et  al.  (1960)  extended 
the  conventional  equation  to  the  form: 

o'  -  (o-u.)  *  .  (2.17) 

Hera  u,  and  u«  represent  pore  air  and  pore  water  pressures  respectively, 
and  %  is  an  empirical  parameter  representing  the  proportion  of  the  soil 
suction  (u.-u«)  that  contributes  to  the  effective  stress.  For  fully 
saturated  soils  x  i*  unity,  and  for  dry  soils  x  zero.  Bishop  et  al. 
(1960)  suggests  that  the  value  of  x  depends  mainly  on  the  degree  ox 
saturation  and  to  a  lesser  degree  soil  structure,  cycle  of  wetting  and 
drying,  and  external  stress  changes.  Blight  (1967)  stated  that  the  main 
drawback  of  Equation  (2.20)  is  the  difficulty  of  evaluating  x> 

3.  Influence  of  Capillarity 

Knowledge  of  the  configuration  of  air  and  water  in  the 
interstices  of  the  soil  mineral  skeleton  helps  to  understand  the  influence 
o*  cupillarity  on  effective  stiress.  Height,  of  capillary  rise  in  soil 
above  the  water  table  dupends  mainly  on  effective  grain  diameter,  whether 
the  water  is  draining  or  imbibing,  and  to  some  extent  particle  shape 
(Holtz  and  Kovacs,  1981).  Pcre  water  above  the  groundwater  table  in 
unsaturated  soil  has  a  negative  pressure  with  respect  to  air  pressure. 
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Bulking  of  •oil  is  a  result  of  capillarity  and  forms  a  very  loose 
relative  density  soil  structure.  Holtz  and  Rovaca  (1981)  describe  how 
sands  under  certain  deposition  conditions  can  form  a  honeycomb  structure 
as  a  result  of  bulking  (Figure  2.8).  Individual  soil  particles  in  this 
unsaturated  soil  structure  are  held  together  by  capillary  stress.  The 
capillary  stress  is  formed  at  the  air  and  water  interface  and  results  in 
interparticle  stress  (Figure  2.9). 

McWhorter  and  Sunada  (1977)  describe  capillary  stress,  or 
interfacial  tension  (f)  between  soil  particles  as  a  force  that  opposes 
pressure  differences  in  the  pore  water  and  local  air  pressure. 
Interfacial  tension  is  measured  as  force  per  unit  length  acting  along  the 
perimeter  of  the  interface  in  a  direction  tangent  to  the  curved  water 
surface. 

Metric  suction,  or  capillary  pressure  p,  as  defined  by  McWhorter  and 
Sunada  (1977)  is  the  difference  between  the  air  and  water  pressure  (u.-u^) 
in  the  soil  pore  spaces.  Capillary  pressure  and  interfacial  tension  are 
related  by 


where  r«  represents  the  radius  of  the  curved  water  interface  (Figure  2.9). 
Equation  2.18  predicts  that  the  capillary  pressure  increases  with 
decreasing  saturation  because  of  a  reduction  in  r„. 

4.  Shear  Strength  and  Volume  Change  of  Unsaturated  Soil 

The  science  for  unsaturated  soils  was  slow  to  develop  until 
its  appearance  by  Bishop  in  the  late  1950 's.  Fredlund  (1985)  sunsnarizes 
shear  strength  and  volume  change  eqpiations  of  unsaturated  soil. 

Terzaghi's  (1943),  classic  shear  strength  equation  for  saturated 
soil  is  written  in  terms  of  effective  stress  parameters  from  the  effective 
Mohr-Coulomb  failure  envelope 
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T  ■  c'  +  (o„-uJ  tan  4' 


(2.19) 


where 

r  "  shear  stress; 

c'  ■  effective  cohesion  intercept; 

o,  •  total  normal  stress; 

u,  >  pore  water  pressure;  and 

■  affective  angle  of  internal  friction. 

Fredlund  (1981)  revised  the  shear  strength  equation  for  unsaturated  soil 
to  the  form 

X  •>  c'  *  {u,-ujt&n  ♦  (o„-u,)  tan  <>'.  (2.20) 

where  u.  is  the  pore  air  pressure,  and  41*  is  the  angle  of  shear  strength 
increase  with  an  increase  in  matric  suction,  (u.-u«,).  Values  for  are 
consietently  less  than  and  are  on  the  order  of  IS  degrees  (Fredlund, 
198S) . 

Equation  (2.20)  can  be  interpreted  as  having  two  cohesion  terms. 
The  second  of  these  cohesion  terms,  apparent  cohesion,  is  a  function  of 
metric  auction,  (u,*-u.)tan  Apparent  cohesion  is  represented  by  a  third 
dimensicnal  extension  of  the  Mohr-Coulorob  failure  criteria  shown  in  Figure 
2.10.  When  matric  suction  goes  to  tero,  apparent  cohesion  becomes  zero 
and  the  shear  strength  is  reverted  back  to  the  saturated  state  of  Equation 
2.19.  Equation  (2.20)  allows  a  smooth  transition  between  saturated  and 
unsaturated  soils  (Fredlund,  1985). 

The  three  dimensional  failure  region  developed  by  Fredlund  (1981)  is 
planar  due  to  a  linear  increase  in  0*’  and  (Figure  2.10).  The  shear 
strength  of  a  particular  soil  can  be  easily  located  on  this  plane  if  (u.- 
u«)  and  (0-u,)  are  knovm. 

The  classic  volume  change  constitutive  relation  in  terms  of 
compressibility  for  saturated  soils  is 
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.10  Thraa^'diaanaional  axtandad  Mohr-Coulomb  failur*  turfaca  (From 
Fradlund,  1985). 
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t  “  Biy  d{a-uj 


(2.21) 


where  e  is  volumetric  strain  «  e.) ,  and  m,  is  the  coefficient  of 
volume  compressibility.  The  coefficient  of  volume  compressibility  is 
simply  the  inverse  of  the  constrained  modulus f  M. 

The  compressibility  constitutive  equation  for  unsaturated  soil  is 

c  ^  Bh'  d(a-Ug)  *  BI2'  d(u^-uj  (2.22) 

where  m,*  is  the  compressibility  of  the  soil  structure  with  respect  to  a 
change  in  (o-u,),  and  m,*  is  the  compressibility  with  respect  to  a  change 
in  (u,-u^)  (Fredlund/  1985). 

Equation  (2.22)  utilizes  matric  suction  similar  to  the  shear 
strength  equation.  The  stress  variable,  again,  y '.elds  a  smooth  transition 
to  the  saturated  case  of  Equation  (2.21)  when  the  absolute  pore  pressure 
is  zero. 

5.  Unsaturated  Soil  Placement 

Holtz  and  Kovacs  (1981)  describe  an  "end-product"  earthwor)c 
specification  that  states  as  long  as  the  contractor  is  able  to  obtain  the 
specified  relative  compaction,  it  doesn't  matter  how  it  is  obtained,  nor 
the  equipment  he  uses.  Earthwork  specifications  for  compaction  of  sands 
are  generally  of  this  type  in  that  only  the  final  dry  density  is 
specified.  Laboratory  teste  to  determine  the  behavior  of  sand  backfill 
are  generally  conducted  in  this  manor  as  well.  The  design  assumption  is 
that  moisture  content  during  compaction  does  not  effect  the  dynamic 
behavior  of  unsaturated  sands. 

C.  CENTRIFUGE  MODELING 

1.  Historical  Background 

The  earliest  engineering  use  of  a  centrifuge  was  introduced  in 
1869,  by  a  Frenchman  named  Phillips,  who  simulated  self-weight  stresses  of 
structural  beams.  His  tests  were  relatively  insignificant,  however,  and 
the  concept  of  increased  gravity  for  modeling  soil  and  rock  did  not  come 
to  its  fruition  until  1931  by  an  American  named  Bucky.  Bucky  developed 
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the  technique  of  eimulating  eelf-welqht  stresses  in  mines.  At  the  same 
time  in  Russia,  Pokrovsky  developed  the  centrifuge  technique  to  determine 
the  stability  of  slopes  in  river  banks.  Through  the  l950's,  centrifuge 
testing  in  the  U.S.  remained  confined  to  mining  applications.  Xn  1966, 
Schofield  and  his  colleagues  at  the  University  of  Cambridge  built  a 
prototype  geotechnical  centrifuge,  in  1969,  Schofield  built  a  1.5  meter 
geotechnical  centrifuge  at  the  University  of  Manchester  Institute  of 
science  and  Technology.  Then  in  1976  Schmidt  and  Schofield  began  dynamic 
work  on  explosive  cratering.  Schmidt  and  Holsapple  are  predecessors  for 
the  development  of  explosive  testing  to  simulate  nuclear  explosives. 

At  this  time  several  centrifuge  facilities  were  located  throughout 
the  U.S.  and  centrifuge  techniques  were  becoming  more  accepted.  Japan, 
Denmark,  Sweden,  Netherlands,  and  Prance  developed  centrifuge  modeling 
facilities  as  well.  In  1979,  modification  of  a  large  centrifuge  began  at 
the  NASA  Ames  Research  Center  in  order  to  make  it  the  largest  capacity 
centrifuge  in  the  U.S. 

Today,  state  of  the  art  centrifuge  testing  is  practical  and  well 
accepted  world  wide  for  most  static  applications.  However,  there  is  a 
considerable  amount  of  research  to  be  conducted  to  further  develop 
simulation  of  dynamic  events. 

2.  Principles  of  Centrifuge  Modeling 

Gravity  effects  are  important  during  scaled  model  testing  in 
order  to  simulate  conditions  that  properly  replicate  the  prototype.  An 
artificial  gravitational  field  is  the  answer  to  this  problem,  and  the 
centrifuge  is  the  most  convenient  tool  to  achieve  this  requirement. 
Prototype  throughout  this  report  refers  to  a  system  that  exists  at  one  g. 

In  order  to  properly  simulate  gravity  induced  stresses,  it  is 
necessary  to  test  an  Nth  scale  model  in  a  gravity  field  R  times  stronger 
(N  g's)  than  that  experienced  by  the  prototype  at  one  g.  Thus,  a  model 
can  be  constructed  of  height  1/N  the  height  of  the  prototype  and  obtain 
similar  overburden  stresses  when  the  model  is  accelerated  at  N  g's.  To 
achieve  required  accelerations  from  the  centrifuge,  the  scaled  model 
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tr«v«lB  in  a  circular  path  with  uniform  circular  velocity,  v,  aa  ahown  in 
Figure  2.11.  Uniform  circular  velocity  of  a  point  in  the  model  over  a 
period  of  time,  At,  is  given  by 


V  = 


e  f 

At 


(2.23) 


where  6  ia  the  angle  displaced  over  At,  and  r  ia  the  radiua  of  the  path  of 
a  point  in  the  model.  Pointa  in  the  model  are  accelerated  in  a  direction 
toward  the  axis  of  the  centrifuge  with  the  magnitude 


a 


c 


r 


(2.24) 


where  a^  ia  the  centripetal  acceleration.  Relative  to  the  model  itself  is 
an  acceleration  acting  away  from  the  axis  equal  in  magnitude  and  ia  termed 
centrifugal  acceleration.  The  number  of  g’a  (N)  at  which  the  model  ia 
accelerated  can  be  expreaaed  ae 


N  = 


(2.25) 


a.  Buc)cingham  Pi  Theory 

The  phyaical  pararoetera  of  a  model  at  a  particular  g 
level  and  the  reapective  parameters  of  a  prototype  are  related  through 
principles  luiown  aa  scaling  laws.  Scaling  laws  are  necessary  sc  that 
model  parameters  can  be  extrapolated  to  represent  prototype  performance 
(Kline,  1965).  These  scaling  laws  are  developed  through  the  principles  of 
dimensional  analysis  and  utilise  a  method  )inown  aa  the  Buclcingham  Pi 
Theory  (Buckingham,  1915). 

The  Buckingham  Pi  Theory  is  a  powerful  technique  that 
non-dimensional izes  the  parameters  of  a  phyaical  system.  These 
non-dimensional  quantitiea  are  known  as  pi  {n)  terms  and  can  be  related 
outside  the  context  of  the  physical  system  with  similar  quantities  of 
another  physical  system  (i.e.  relating  the  physical  system  of  a  model  with 
that  of  a  prototype).  Scaling  laws  are  then  developed  which  relate  the 
two  systems,  scaling  relations  for  centrifuge  modeling  are  given  in  Table 
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Figure  2.11  Circular  path  of  uniform  circular  velocity  on  the  cen 


2.2.  Derivation  of  scaling  laws  and  further  detail  of  the  Buckingham  Pi 
Theory  are  provided  in  APPENDIX  A. 

b.  Modeling  of  Models 

Scaling  laws  can  become  very  complicated  when  several 
parameters  interact  with  one  another  resulting  in  distortion  of  the  model. 
In  order  to  prevent  this,  one  must  verify  the  scaling  assumptions,  and 
this  is  best  accomplished  by  full  scale  test  .ng  of  the  prototype. 
Prototype  testing  can  be  very  costly,  however,  and  sometimes  impossible. 
As  an  alternative,  the  concept  of  "modeling  of  skodels"  provides  a  check  on 
the  consistency  of  centrifuge  modeling. 

Modeling  of  models  is  illustrated  by  Ko  (1988)  in  Figure  2.12  where 
the  model  size  is  plotted  against  the  gravity  level  on  log-log  scale. 
Consider  a  1000  cm  prototype  at  point  Al.  If  the  prototype  is  modeled  at 
1/10  scale,  the  dimensions  become  100  cm  at  10  g's  (A2),  or  at  1/100  scale 
the  dimensions  become  10  cm  at  100  g's  (A3).  Points  A2  and  A3  are  models 
of  Al  and  are  also  models  of  each  other.  To  verify  that  the  performance 
of  the  models  simulates  the  prototype,  A2  and  A3  can  be  compared  to  one 
another  in  the  absence  of  the  prototype. 

3.  Limitations  of  Models 

Before  conducting  model  testing,  the  applicable  scale  to  be 
used  must  be  evaluated  outside  of  the  scaling  laws  and  theory  of  multiple 
gravity  tasting.  The  physical  limitations  of  the  modeling  material  plays 
a  critical  role  in  the  determination  of  scale.  Using  the  same  soil 
material  to  construct  both  the  model  and  prototype  is  a  great  advantage 
when  conducting  centrifuge  modeling.  By  doing  this,  the  complex  modeling 
of  constitutive  soil  properties  can  be  avoided.  Furthermore,  homologous 
points  in  the  geometrically  similar  model  and  prototype  will  be  subjected 
to  similar  stresses,  and  thus  will  develop  the  same  strains  (Schmidt, 
1981;  and  Ko,  1988). 
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Tabla  2.2  STANDMU)  SCALING  LAWS  (FROM  BRADLEY  BT  AL. ,  1984). 


N  -  NUMBER  OF  G  ACCBLERATIONS 


QUANTITY  PROTOTYPE 

LINEAR  DIMENSION  1 

GRAVITY  (g)  1 

AREA  1 

VOLUME  1 

DYNAMIC  TIME  1 

VELOCITY  (DISTANCE/TIME)  1 

ACCELERATION  (DISTANCE/TIME^)  1 

DENSITY  (MASS /VOLUME)  1 

UNIT  WEIGHT  (FORCE/UNIT  VOLUME)  1 

FORCE  1 

STRESS  (FORCE/ AREA)  1 

MASS  1 

ENERGY  1 

STRAIN  (DISPLACEMENT/UNIT  LENGTH)  1 

HYDRODYNAMIC  TIME  1 

IMPULSE  1 


MODEL 

1/N 

N 

1/N’ 

1/N’ 

1/N 

1 

N 

1 

N 

1/N’ 

1 

1/N’ 

1/N’ 

1 

1/N’ 

1/N’ 
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Sixa  of  Kodol 


1  10  100  1000 


Flguxo  2.12  Or»phic»l  s«pr»»«nt*tion  of  ‘Bodoliiag  of  nodoli*  (Xo,  1988). 
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Although  scaling  theory  suggests  that  particle  site  should  be  scaled 
with  g's,  perturbations  from  similar  model-prototype  materials  have  proven 
insignificant  as  long  as  the  grain  site  is  maintained  much  smaller  than 
adjacent  structure  elements.  Bradley  et  al.  (1984)  reports  that  a  siaiple 
reduction  in  particle  site  has  proven  unsatisfactory  due  to  the  complex 
interrelationship  of  soil  strength  and  response  to  particle  site.  Testing 
based  on  "modeling  of  models"  can  provide  guidelines  for  limitations 
arising  from  grain  site  effects. 

Model  constructability  depends  upon  the  limitations  of  the 
centrifuge  facility.  As  discussed  earlier,  centrifugal  acceleration  in 
the  model  is  dependent  on  radial  distance  from  the  center  of  rotation. 
Thus,  stresses  vary  nonlinearly  with  depth  in  the  model  because  of  the 
difference  in  centrifugal  force  with  depth.  This  can  complicate  model 
behavior  unless  the  radius  arm  to  the  model  is  sufficiently  large  with 
respect  to  model  depth.  A  similar  problem  exists  with  a  difference  in 
stresses  in  a  horizontal  plane  of  the  model  due  to  differences  in  radial 
distance  from  the  center  of  rotation.  Again,  stress  differences  are 
minimised  with  a  sufficiently  large  centrifuge  radius. 

Finally,  scaling  of  instrumentation  sire  and  mass  is  just  as 
important  as  scaling  of  the  model  itself.  Instruments  used  to  measure 
modeling  events  are  often  too  large  to  adequately  represent  prototype 
behavior  at  homologous  points.  Therefore,  it  is  imperative  that  the 
presence  of  these  instruments  do  not  adversely  affect  the  model  behavior. 
Miniature  transducers  are  needed  to  obtain  a  soil-instrument  likeness  in 
scaled  size. 
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SECTION  III 


IXPERIKENTAI.  INVESTIGATION 

A.  CENTRIFUGE  FACILITY,  INSTRUMENTATION  AND  EXPLOSIVES 

1.  Dttscription  of  the  C«ntrifuga 

The  centrifuge  used  to  conduct  the  experimental  phase  of  this 
research  is  a  Genisco  (model  E185,  Serial  Number  11)  located  at  HQ 

AFCESA/RACS,  Tyndall  Air  Force  Base,  Florida  (Figure  3.1).  The  capacity 

of  the  centrifuge  is  13.6  metric  g-tons  at  a  maximum  acceleration  of  100 

g's  and  radius  of  1.83  meters.  Thus,  at  its  maximum  acceleration  of  100 

g's,  a  load  of  136  kg  can  be  applied  to  each  of  the  two  payload  platforms. 
The  centrifuge  is  rotated  by  a  variable  speed  hydraulic  motor  regulated 
from  the  control  console  shown  in  Figurs  3.2. 

The  centrifuge  consists  of  two  symmetrical  cantilever  arms  opposite 
of  one  another,  each  of  which  supports  a  76  cm  square  payload  platform. 
The  specimen  payload  is  attached  to  one  platform  while  counter  balance 
weight  of  eqpjal  mass  is  placed  on  the  opposite  platform.  While  the 
centrifuge  is  in  flight,  an  automatic  dynamic  balancing  motor  equalizes 
payload  masses  of  loss  than  4.53  kg.  The  rotation  rate  of  the  centrifuge 
is  indicated  on  a  digital  tachometer  locatsd  on  the  control  console 
(Figure  3.2). 

2.  Instrumentation 

a.  Carbon  piezoresistive  stress  gages 

Stress  time  histories  generated  by  the  blast  are 
measured  with  1/6  watt  1000  ohm  {*S%  tolerance)  carbon  resistors 
manufactured  by  Allen  Bradley  (Figure  3.3).  These  resistors  were  selected 
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Figure  3.1  Centrifuge  located  at  Tyndall  Air  Force  Base. 
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for  their  small  size  (1.58  mm  in  diameter  and  4.00  on  in  length)  and  their 
ability  to  represent  a  point  measurement.  In  addition,  carbon  resistors 
are  readily  available,  can  be  statically  calibrated,  and  are  inexpensive. 

Carbon  stress  gages  are  configured  as  Wheatstone  bridges  and  have 
four  carbon  resistors  (Figure  3.4).  The  active  arm  of  the  gage,  R,,  is 
located  in  the  specimen,  while  the  remaining  bridge  arms  are  mounted 
outside  the  specimen  on  the  centrifuge.  All  resistors  are  nominally  equal 
to  one  another  (i.e.  R|«R]>R)*R,) .  Operation  of  the  bridge  utilizes  the 
stress-resistance  characteristics  of  carbon  such  that  a  resistance  change 
due  to  applied  stress  in  the  active  arm  produces  a  proportional  voltage 
change  in  the  bridge. 

As  described  by  Holloway  et  al.  (1985),  carbon  resistors  display  a 
very  fast  response  time,  roughly  the  time  required  for  the  shocli  wave 
front  to  cross  the  resistor.  Carbon  resistors  can  be  used  to  measure 
dynamic  peak  pressures  to  at  least  100  MPa. 

Based  on  the  agreement  between  static  and  dynamic  calibration 
(Holloway,  1985),  quasi-static  stress-voltage  relationships  for  eight 
carbon  stress  gages  were  experimentally  determined.  Linear  stress-voltage 
relationships  were  obtained  for  these  gages,  and  a  slope  of  the  mean 
regression  of  190.7  MPa/volt  resulted.  All  voltage-time  history  data 
analyzed  in  this  report  utilizes  this  calibration  value.  Gage  calibration 
results  are  given  in  appendix  B. 

b.  Accelerometers 

The  instruments  used  to  determine  peak  free  field 
accelerations  are  BNDEVCO  piezoresistive  accelerometers  (model  7270A) 
shown  in  Figure  3.5.  These  accelerometers  are  rugged  undunped  units 
designed  for  linear  shock  measurements  up  zo  20,000  g's.  The  units  are 
chemically  sculptured  from  a  single  piece  of  silicon  with  an  active  four 
arm  strain  gage  Whsatstons  bridge.  Its  low  mass  of  l.S  grams,  small  size 
(14.2  mm  by  7.1  rom  by  2.8  am  thick),  high  resonant  frequency  (mounted;  350 
kHz),  and  zero  damping  allow  the  acceleromatars  to  respond  accurately  to 


35 


Lgure  3.5  Endevco  Piezoresi-stivs  accelerometers . 
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fast  rite  time,  short  duration  shock  motion.  Calibration  of  the 
accelsrometera  was  conducted  by  ENDEVCO.  Calibration  curves  are  given  in 
APPENDIX  B. 


c.  Data  Recording  System 

Voltage-time  histories  %Mre  recorded  by  a  Portable  Data 
Acquisition  System  (Model  S700)  manufactured  by  Pacific  Instruments.  It 
conditions,  amplifies,  digitises,  and  records  transducer  signals  at 
sampling  rates  of  up  to  one  million  samples  per  second.  The  fully 
programmable  16  channel  system  has  a  memory  capacity  of  256  kilobytes  per 
channel.  The  data  acquisition  system  is  mounted  on  the  axis  of  the 
centrifuge  (Figure  3.1)  so  that  data  can  be  transferred  directly  from  the 
instruments  without  having  to  send  signals  through  the  slip  rings.  Slip 
rings,  located  at  the  base  of  the  centrifuge  axis,  provide  all  electrical 
communication  to  the  control  console  during  operation.  A  high-speed  video 
camera  is  mounted  above  the  data  recorder  (Figure  3.1)  and  allows  the 
operator  to  continuously  monitor  the  specimen  during  flight. 

3.  System  Configuration 

A  schematic  of  the  centrifuge  and  equipment  is  shown  in  Figure 
3.6.  When  the  firing  mechanism  is  engageo,  a  signal  is  sent  through  the 
slip  rings  triggering  the  data  acquisition  and  detonating  the  charge 
simultaneously.  After  tecting  is  completed,  an  IEEE-488  general  purpose 
interface  bus  (GPIB)  ie  connected  to  the  data  recorder  mounted  on  the 
centrifuge,  and  the  data  is  downloaded  to  a  personal  computer.  A  system 
language,  ASYST,  is  utilised  for  communication  between  the  data  recorder 
and  computer  by  means  of  the  GPIB.  Software  utilised  to  control  data 
recorder  functions  and  channel  progranvning  is  Pacroon  revision  C  written 
by  Pacific  Instruments.  The  large  amount  of  data  acquired  over  a  period 
of  very  fast  sampling  is  stored  on  a  44-roegabyte  removable  cartridge 
Bernoulli. 
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Centrifuge 

Model 


rigur*  3.6  Schematic  of  th«  contrifuga  *nd  aquipraant. 
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4.  Explosivtts 

The  explosive  charges  utilised  in  these  tests  are  exploding 
bridgewire  detonators  manufactured  by  Reynolds  Industries.  Two  types  of 
detonators  were  used  in  our  studies,  the  KP~H2  detonator  and  a 
modification  of  this  detonator,  h  cross  section  of  an  RP-83  detonator  is 
shown  in  Figure  3.7. 

The  detonators  consist  of  an  exploding  wire  bridge,  80  mg  of  a  low- 
density  initiating  explosive  of  PETN  (pentaerythrito?  tetranitrate) ,  123 
mg  of  PBX  9407  (cyclotetramethyl-enetetranitramine)  in  addition  to  the 
initiator,  and  a  high  density  PBX  9407  output  charge,  all  of  which  are 
contained  in  a  0.018  mm  thick  aluminum  cup.  The  high  output  charge  of  the 
RP-83  detonator  consists  of  four  individual  preesings  of  PBX  9407  each 
weighing  227  mg,  and  the  modified  RP-83  detonator  has  only  one  pressing 
of  PBX  9407.  Table  3.1  breaks  down  the  total  output  charge  for  both 
detonators. 

The  composition  of  PBX  9407  consists  of  94  percent  ROX,  and  6 
percent  Exon  461  which  acts  as  a  binder.  Explosive  TNT  equivalency  for 
ROX,  PETN  and  PBX  9407,  as  described  in  The  Manual  for  the  Prediction  of 
Blast  and  Fragment  loadings  on  Structures  (Baker  et  al,  1980),  is  given 
in  Table  3.2.  Most  blast  work  is  presented  using  ecguivalent  TNT  weight. 
TNT  equivalents  for  shock  presoures  are  based  on  heat  of  detonation  as 
suggested  in  Anonymous  (1990). 

Because  the  charge  density  of  the  PETN  initiator  is  about  half  of 
that  of  PBX  9407  (Table  3.1),  and  because  its  mass  is  small  compared  to 
the  total  output  charge,  the  mass  of  PETN  is  not  included  in  the  total 
output  charge  of  the  detonator.  These  total  output  charge  masses  were 
also  suggested  by  Reynolds  Industries  (personal  communication,  Ron 
Varosh).  Future  reference  to  detonator  explosive  mass  in  this  report  will 
refer  to  PBX  9407  total  output  charge  (Table  3.1^. 
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PARTS  DESCRIPTION 

1.  MOLDED  HEAD; 

2.  STRING. 

3.  BRIOGEWIRE;  GoW. 

4.  INmATING  EXPLOSIVE;  80  mg  Of  PETN. 

5.  TOTAL  OLTTPUT  CHARGE;  1031  mg 

6.  ALUMINUM  CUP: 


?igur«  3.7  RP-83  datonator  (Aaynolda  Induatriaa). 


Table  3.1  CHARGE  MASS  FOR  RP-83  AND  MODIFIED  RP-83  DETONATORS 

(PERSONAL  COMMUNICATION:  RON  VAROSH,  REYNOLDS  INDUSTRIES). 


CHARGE  DETONATOR 


COMPONENT 

TYPE 

DENSITY 

RP-83  Modified  RP-83 

(gm/cc) 

(mg) 

(mg) 

Initiating 

Exploaive 

PETN 

0.88 

80* 

80* 

Pallet 

PBX  9407 

1.60 

123 

123 

Output  Charge 

PBX  9407 

1.60 

908 

227 

TOTAL  OUTPUT 

CHARGE: 

1031  mg 

350  mg 

*  Not  uaed  in  total  output  charge. 


Table  3.2  EXPLOSIVE  EQUIVALENCE  WITH  RESPECT  TO  HEAT  OF  DETONATION  (BAKER 


ET  AL. ,  1980) . 

HEAT  OF 

DETONATION 

TNT 

EQUIVALENCE 

EXPLOSIVE 

(ft-lb/lb) 

TNT 

1.97xl0‘ 

1.000 

PETN 

2.31x10' 

1.173 

RDX 

2.27x10* 

1.152 

PBX  9407 

2.24x10* 

1.137 
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B.  SOIL  PROPERTIES  RND  SPECIMEN  PREPARATION 

1.  Soil  Propart 

Tyndall  baach  aand  ia  obtainad  from  dunaa  at  Tyndall  Air 
Forca  Baaa,  Florida.  Tha  aand  ia  ovan  driad,  than  aiavad  to  ramova 
organica.  Ottawa  20-30  aand,  ia  obtainad  from  tha  U.S.  Silica  Company, 
Ottawa,  lllinoia  (A8TH  C190).  A  auaaaary  of  tha  phyaical  propartiaa  of 
both  Tyndall  and  Ottawa  aanda  ia  givan  in  Tablaa  3.3a  and  3.3b.  Grain 
aiza  diatributiona  for  tha  aanda  are  ahown  in  Figura  3.8  (ASTM  D421-56  and 
D422-63).  Both  Tyndall  and  Ottawa  aanda  ara  poorly  gradad  aanda. 

Deaaturation  curvea  obtainad  by  the  poroua  plate  method  for  Tyndall 
and  Ottawa  aanda  ara  ahown  in  Figuraa  3.9(a)  and  (b)  (ASTM  D2325).  Both 
aanda  begin  to  daaaturata  at  capillary  praaauraa  (u.-u..)  of  1  ]cPa  and 
reach  raaidual  aaturationa  (approximately  11  percent  for  Tyndall  aand  and 
5  percent  for  Ottawa  aand)  at  capillary  praaauraa  graatar  than  15  kPa. 
ASTM  atandard  teat  mathoda  can  be  found  in  ASTM  ( 1987 ) . 

2.  Compaction  Mathoda 

Soil  placement  la  critical  in  conducting  our  invaatigation  and 
wa  ara  intereatad  in  placement  tachniquaa  aimilar  to  thoaa  conntonly  uaad 
to  place  coheaionlaaa  backfill  material  during  conetruction.  Two  methods 
of  compaction,  raining  and  vibration,  were  avaluatad  for  preparing  aand 
apaclroena  efficiently  and  with  reproducibility.  The  vibration  method  waa 
uaad  to  conduct  thia  caaeareh  effort. 

The  method  of  raining  ia  widely  uaad  and  offara  tha  ability  to 
obtain  aand  apecimena  with  a  great  deal  of  reproducibility  and  beat 
simulataa  tha  natural  procaaa  of  aand  depoaition.  Tha  procedure  for  thia 
method  conaiata  of  dry  aand  raining  from  a  ahuttar,  through  a  vertical 
column,  than  through  a  aeries  of  diffuaar  siavaa.  Shutter  porosity  and 
fall  height  ara  a  few  of  tha  many  variables  one  can  adjust  to  obtain  a 
target  density.  Three  apparent  problems  exist  with  this  mathodi  first, 
there  are  no  )cnown  techniques  for  raining  partially  saturated  apecimena; 
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Table  3.3a  PHYSICAL  PROPERTIES 

OF  TYNDALL  BEACH  SAND. 

Soil  Daacription 

Poorly  graded  aand  compoaed 
mainly  of  quartz.  Soil 
particlea  are  aubangular  to 
aubrounded. 

Maxinuffl  Relative  Denaity 

1630  kg/m’ 

(ASTM  D4253) 

Minimum  Relative  Denaity 

1450  kg/m’ 

(A8TM  D4254) 

Maan  Grain  Size  (Ojo) 

0. 25  nm 

(ASTM  D422) 

Specific  Gravity  (G,) 

2.65 

(ASTM  D854) 

Table  3.3b  PHYSICAL  PROPERTIES 

OF  OTTAWA  20-30  SAND. 

Soil  Description 

Poorly  graded  sand  composed 
mainly  of  quarts.  Soil 
particles  are  aubrounded 
to  rounded  in  shape. 

Maximum  Relative  Density 
<P4..U».) 

1720  kg/ro’ 

(ASTM  D4253) 

Minimum  Ralative  Density 

(PdiriiiH.) 

1560  kg/m’ 

(ASTM  D42S4) 

Mean  Grain  Size  (Djg) 

0.70  mm 

(ASTM  D422) 

Specific  Gravity  (0,) 

2.65 

(ASTM  0654) 
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PERCENT  PASSING  BY  WEIGHT 


Figurti 


3 >8  Crain  aiz*  discribucions  for  Tyndall  baacb  and  Ottawa  20-30 

aands. 
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CAPILLARY  PRESSURE  (kPa) 


■•cond,  placwnent  of  a  level  layer  ie  difficult  to  obtain;  and  third, 
raining  is  not  used  for  placing  backfill  during  construction. 

Vibration  proves  to  be  quite  effective  for  compacting  both  dry  and 
wet  soils  and  opposing  the  forces  of  capillarity.  The  vibration  apparatus 
used  for  preparing  the  specimens  is  a  variable  amplitude  60  He  vibrator 
oftotor,  manufactured  by  Syntron  (model  V51  Dl),  and  is  mounted  to  a  13  mm 
thick  aluminum  plate  (figure  3.10).  The  apparatus  is  placed  on  the 
surface  of  the  soil  and  allowed  to  compact  as  a  vibrating  surcharge. 
Additional  surcharge  weight  can  be  placed  on  the  aluminum  plate.  The 
vibration  technique  developed  is  similar  to  vibratory  methods  used  during 
construction  (i.e.  walk  behind  vibrators  and  vibrating  roller  compactors). 

3.  Specimen  and  Instrumentation  Placement 

Centrifuge  soil  models  are  constructed  in  a  specimen  bucket 
with  inside  diameter  of  46  cm  and  depth  of  25.7  to  27.9  cm  (figure  3.11). 
Models  consist  of  five  horizontal  lifts,  each  lift  is  compacted  separately 
in  order  to  obtain  uniformity  throughout  the  specimen.  Target  degree  of 
saturation  for  the  specimens  are  obtained  when  the  appropriate  amounts  of 
sand  and  water  for  each  lift  are  combined  and  compacted  to  a  unit  lift 
volume. 

Target  saturation  levels  for  Tyndall  sand  are  0,  17,  35,  53,  and  70 
percent,  and  for  Ottawa  sand  are  0,  20,  40,  and  60  percent.  Target  dry 
density  for  Tyndall  sand  is  1521  kg/m’  (95.0  pcf)  and  for  Ottawa  sand  is 
1612  kg/m’  (100.7  pcf).  Dry  relative  densities  for  Tyndall  and  Ottawa 
sands  ars  42  percent  and  35  percent  respectively.  Distilled  water  is  used 
as  the  pore  fluid. 

The  bottom  four  lifts  of  the  specimen  are  5.1  cm  thick  and  the 
overburden  lift,  above  the  instrument  layer,  varies  in  height  (figure 
3.11).  Overburden  stresses  at  a  prototype  depth  of  1.43  meters  are 
simulated  with  a  model  overburden  of  7.6  cm  when  accelerated  to  16.86  g's. 
Similarly,  accelerating  the  model  to  26.34  g's  with  an  overburden  of  5.4 
cm  will  produce  a  eimilar  prototype  overburden  etrees.  As  noted  in 
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Overburden 

Lift 


Lift  4 


u';.  3 


Lift  2 


Lift  1 

K  ■  -  -■  ■" 


Charge 

4- 


Soedmen 
Bucket 


Instrumentation 

— x X ^ 


.s 

N 

V 

V 
■V 

V 

V 


_i_ 

5.4-7.6  cm 


20.3  cm 
(5.icm/lift) 


t 


46  cm 


Figure  3.11  Croee  sect ion  of  centrifuge  epecimen. 
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Section  II. A. 3,  th*  depth  of  burial  of  the  blast  is  critical  in  order  to 
completely  transmit  the  energy  of  the  blast  to  the  surrounding  medium. 
A  scaled  depth  of  burst  of  0.71  m/kg*'^  (1.8  ft/lb'^)  was  chosen  to  insure 
a  ground  shock  coupling  factor*  t,  of  unity  (Figure  2.1).  Tables  3.4a  and 
3.4b  list  model-prototype  relations  for  Tyndall  and  Ottawa  sands. 
Kxtrapolating  the  prototype  performance  using  scaled  models  at  different 
g  levels  allows  one  to  verify  the  principle  of  "modeling  of  models." 

At  the  instrument  layer,  carbon  stress  gages  and  accelerometers 
extend  outward  form  the  center  of  the  bucket  where  the  charge  is  placed. 
Figure  3.12  illustrates  the  gage  positions  and  distances  from  the  charge. 
The  prefix  letter  R  designates  carbon  stress  gages  and  the  letter  A 
designates  accelerometers.  Gage  wires  trail  behind  the  gages  toward  the 
bucket  wall  so  that  stress  waves  traveling  to  the  gages  will  not  be 
distorted.  Data  interference  due  to  transient  waves  reflected  off  the 
bucket  bottom  or  transmitted  through  the  container  are  minimized  by 
situating  the  gages  in  the  path  of  shortest  travel  from  the  blast  (Figure 
3.13). 


4.  Charge  Placement 

Blast  simulation  requires  that  the  expxosive  yield  of  a  model 
is  the  reciprocal  of  the  acceleration  cubed  (1/N’)  according  to  the 
scaling  law  for  mass  in  Table  2.2.  Depth  of  burial  of  the  explosive 
similarly  scales  with  the  reciprocal  of  acceleration  using  the  linear 
dimension  property.  Model-prototype  relationships  for  charge  mass  and 
depth  of  burial  for  Tyndall  and  Ottawa  sands  are  given  in  Tables  3.4a  and 
3.4b. 

Two  methods  are  used  for  placing  charges  consistently  and  with 
minimum  disturbance  to  the  model t  first,  a  drilling  method  for  placement 
in  moist  specimens,  and  second,  a  vacuum  method  for  dry  specimens.  The 
depth  of  the  borehole  is  determined  such  that  the  centroid  of  the  charge 
is  at  the  same  level  as  the  instrument  layer.  Dry  sand  is  backfilled  into 
the  borehole  after  the  charge  is  placed. 


SO 


Table  3.4a  HODEL-PROTOTYPE  CHARACTERISTICS  FOR  TYNDALL  BEACH  SAND 


HODEL  PROTOTYPE 


Amt 

Accel. 

Charge 

Overburden 

Charge 

Overburden 

RPH 

Radius 

N* 

W^ 

D 

N*xW 

NxD 

(m) 

(g'«» 

(gin) 

(cm) 

(kg) 

(m) 

102.1 

1.83 

18.86 

1.031 

7.6 

6.9 

1.43 

120.6 

1.83 

26.34 

0.35 

5.44 

6.4 

1.43 

*  Acceleration 

at  level 

of  instrumentation 

(r  -  1.62 

m) . 

^  Baaed 

on  total  output 

charge 

mass  (Table 

3.1). 

Table  3.4b  MODEL-PROTOTYPE  CHARACTERISTICS  FOR  OTTAWA  20-30  SAND. 


MODEL 

PROTOTYPE 

RPM 

Arm 

Radius 

(m) 

Accel. 

N* 

<g'») 

Charge 

W* 

(gin) 

Overburden 

D 

(cm) 

Charge 

N’xW 

(kg) 

Overburden 

NxD 

(m) 

191.9 

1.83 

66.65 

0.35 

5.44 

104 

3.63 

*  Acceleration  at  level 

*  Based  on  total  output 

of  instrumentation 
charge  mass  (Table 

(r  ■  1.62 
3.1) . 

m) . 
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Tigur*  3.12  Top  viaw  et  gage  and  charge  placement  in  eentri 
■pecimena . 


P'lmtry 
t^tvf  Ptm 


Figure  3.13  Illuetration  of  wave  path  in  centrifuge  apeeiaen. 
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SECTION  IV 


EXPERIMENTAL  RESULTS 

This  chapter  presents  the  results  of  10  centrifuge  explosive  tests 
on  Tyndall  beach  sand  and  4  centrifuge  explosive  tests  on  Ottawa  20-30 
sand. 


A.  STRESS  GAGE  RESULTS 

Voltage-time  histories  for  four  typical  carbon  stress  gages  are 
illustrated  in  Figure  4.1.  Voltage  readings  were  converted  to  values  of 
stress  using  the  calibration  constant  190.7  HPa/volt  obtained  from  static 
calibration  of  eight  carbon  gages  {APPENDIX  B).  The  stress  gages  are 
located  at  increasing  distances  from  the  charge  (Figure  3.13).  Arrival 
time  of  the  stress  wave  increases  respectively  with  increasing  distance 
to  the  gages.  Attenuation  of  pea)(  stress  is  evident  from  the  decrease  in 
wave  amplitudes  with  time  in  Figure  4.1. 

A  summary  of  stress  gage  results  for  Tyndall  Beach  sand  are 
presented  in  Table  4.1  listing  values  of  peak  stress,  wave  arrival  tiroes, 
and  average  wave  velocity  between  respective  stress  gages.  Ten  tests  were 
conducted,  two  tests  each  of  0,  20,  40,  60,  and  70  percent  s<.turation8. 
Test  names  are  listed  in  the  left  column  of  these  tables;  the  first 
numeral  represents  the  g  level  at  the  model  platform,  followed  by  the 
saturation  of  the  test,  then  the  teet  number  performed  at  that  saturation 
level.  The  g  level  at  the  instrumentation  level  is  given  in  parenthesis 
below  the  test  name.  Results  of  four  tests  on  Ottawa  20-30  sand  are  given 
lit  Table  4.2  for  saturations  of  0,  20,  40,  and  60  percent. 

Peak  stresses  from  Tables  4.1  and  4.2  for  Tyndall  and  Ottawa  sands 
are  plotted  versus  scaled  distance  in  Figures  4.2  and  4.3.  Wave  arrival 


(The  reverse  of  this  page  is  blank.) 
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Tabl*  4.1  8TRBSS  OAGB  RKSULTS  FOR  TITNDRLL 


Test 

Charge 

WL 

(gm) 

Instr. 

Gage 

Oist.  to 
Charge 
(cm) 

Peak 

Stress 

(kPa) 

Arrival 

Time 

(us) 

21  3-0-7 

1.031 

R1 

3.81 

17015 

47 

(1b.86) 

R2 

7.62 

3584 

114 

R3 

12.7 

132 

221 

R4 

17.78 

528 

352 

RS 

20.32 

183 

436 

R9 

17.78 

416 

3*' 4 

RIO 

20.32 

214 

483 

29.8-0-9 

0.35 

R1 

3.81 

11172 

47 

(26.34) 

R2 

7.62 

2030 

115 

R3 

12.7 

410 

260 

R4 

17.78 

130 

464 

RS 

20.32 

114 

565 

R6 

17.78 

336 

413 

fl7 

20.32 

86 

532 

29.8-20-2 

0.35 

R1 

3.81 

7464 

51 

(26.34) 

R2 

7.62 

2978 

122 

R3 

12.7 

601 

249 

R4 

17.78 

370 

377 

R5 

20.32 

109 

449 

21.3-20-3 

1.031 

R1 

3.81 

15917 

42 

(18.86) 

R2 

7.62 

4004 

114 

R3 

12.7 

1243 

223 

R4 

20.32 

463 

426 

RS 

17.78 

536 

348 

R9 

17.78 

360 

368 

RIO 

20.32 

286 

416 

29.8-40-1 

0.35 

R1 

3.81 

6291 

57 

(26.34) 

R2 

7.62 

1573 

132 

R3 

12.7 

597 

262 

R4 

17,78 

282 

427 

R5 

20.32 

126 

516 

RIO 

20.32 

183 

478 

21.3-40-2 

1.031 

R1 

3.81 

22752 

41 

(18.86) 

R2 

7.62 

4404 

111 

R3 

1^7 

1285 

214 

R4 

17.78 

757 

367 

R5 

20.32 

377 

445 

R9 

17.78 

702 

359 

RIO 

20.32 

429 

409 

BXACH  SRHD 

Wave 

Velocity 

(m/s) 


S69 

475 

388 

302 

257 

560 

350 

249 

251 

213 


537 

400 

397 

353 

529 

466 

375 

326 

529 


508 

391 

308 

285 


544 

493 

332 

326 

508 
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4.1  (oontinusd) 


STRBSS  OAGB  RBSUI.T8  BOR  lYNDALL  BEACH  SAND 


Charge 

Test  Wt  Instr. 
(gm) 


29.8-60-1  0.35  R1 

(26.34)  R2 

R3 

R4 

R5 

RIO 

21.3-60-2  1.031  R1 

(18.86)  R2 


R3 

R4 

RS 

flS 

Rio 

21.3-70-3  1.031  R2 

(18.86)  R3 

R4 

RS 

R6 

29.8-70-4  0.35  R1 

(26.34)  R2 

R3 

R4 

R5 

R6 


Gage 

Oist.  to 
Charge 
(cm) 

Peak 

Stress 

(kPa) 

Arrival 

Time 

(us) 

3.81 

5977 

60 

7.62 

1098 

131 

12.7 

387 

260 

17.78 

290 

433 

20.32 

132 

597 

20.32 

194 

480 

3.81 

20534 

49 

7.62 

3937 

112 

12.7 

1336 

216 

17.78 

694 

362 

20.32 

368 

440 

7.62 

5008 

113 

12.7 

1401 

226 

7.62 

4720 

113 

1^7 

1584 

225 

17.78 

644 

3S8 

20.32 

370 

432 

20.32 

667 

431 

3.81 

7039 

61 

7.62 

1443 

139 

12.7 

467 

283 

20.32 

120 

543 

20.32 

69 

628 

17.78 

137 

443 

Parentheses  indicate  g  level  at  instrument  depth. 


Wave 

Velocity 

(m/s) 


537 

394 

294 

0 


605 

488 

348 

326 

450 


454 

382 

343 


488 

353 

293 

0 
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TAbla  4.2  STKBSS  0A6S  USULTS  FOR  OTTAWA  20-30  SAND 


Test 

Charge 

Wt. 

(gm) 

Instr. 

Gage 
Dist.  to 
Charge 
(cm) 

Peak 

Stress 

(kPa) 

Arrivai 
Time 
(u  s) 

75.3*0-10 

0.35 

R1 

3.81 

8916 

20 

(66.65) 

R2 

7.62 

503 

102 

R3 

12.7 

349 

279 

R4 

17.78 

103 

468 

R5 

20.32 

19 

652 

R6 

20.32 

29 

634 

75.3-20-6 

035 

R2 

7.62 

1192 

120 

(66.65) 

R3 

127 

448 

273 

R4 

17.78 

92 

487 

R5 

20.32 

59 

546 

75.3-40-5 

0.35 

R1 

3.31 

6293 

47 

(66.65) 

R2 

7.62 

1312 

128 

R3 

127 

391 

269 

R4 

17.78 

162 

647 

R5 

20.32 

23 

576 

TS.3-60-5 

0.35 

R1 

.?■  Rt 

4879 

51 

(66.65) 

R2 

7.62 

1333 

138 

R3 

12.7 

381 

289 

R4 

17.78 

72 

520 

R5 

20.32 

27 

618 

R6 

20.32 

74 

644 

Wave 

Velocity 

(m/s) 


46S 

2B7 

269 

138 

153 


332 

237 

431 


470 

260 

274 

248 


438 

336 

220 

259 

205 


Parentheses  indicate  g  level  at  instrviment  depth. 
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tiiMS  from  Tablaa  4.1  and  4.2  were  obtained  front  otreae-tine  hiatoriea  at 
tha  wave  front,  or  the  point  at  which  the  change  in  streea  with  time  ia 
greateat.  Streaa  wave  velocity  between  gagea  waa  determined  by  taking  the 
difference  in  the  gage  poaitiona  from  the  charge,  then  dividing  by  tha 
difference  in  arrival  timea  to  the  respective  gagea.  The  wave  velocity 
of  the  model  acalea  directly  with  the  wave  velocity  of  the  prototype 
(Table  2.2).  Plota  of  atreaa  wave  velocity  veraua  acaled  dietance  for 
Tyndall  and  Ottawa  aande  are  ahown  in  Pigurea  4.4  and  4.5. 

B.  ACCELEROMETER  RESULTS 

Accelerometera  were  utilized  for  free  field  particle 
acceleration-time  hiatoriea.  Calibration  of  the  gagea  waa  conducted  by 
ENDEVCO  and  programmed  into  the  data  recorder  ao  that  valuea  of 
acceleration  in  unite  of  g'a  could  be  directly  obtained  (1  g  *  9.81  m/a^ 
-  32.2  ft/e^). 

A  typical  acceleration-time  hiatory  curve  ia  illuatrated  in  Figure 
4.6.  Peak  particle  accelerationa  were  obtained  from  the  wave  peaka. 
Particle  velocity-time  hiatoriea  were  obtained  by  integrating 
acceleration-time  hiatoriea  (Figure  4.6).  Similarly,  peak  particle 
velocitiea  were  obtained  from  the  velocity-time  wave  peaka. 

Peak  particle  accelerations  and  peak  particle  velocitiea  for  Tyndall 
and  Ottawa  aands  are  summarized  in  Tables  4.3  and  4.4.  Scaled  peak 
particle  accelerations  are  plotted  veraua  scaled  dietance  in  Figures  4.7 
and  4.8.  Peak  particle  accelerations  are  scaled  by  multiplying  by  N  so 
that  prototype  peak  particle  accelerationa  are  simulated  (Table  2.2). 

Peak  particle  velocities  obtained  from  velocity-time  hiatoriea  were 
obtained  in  units  of  g-seconda,  then  converted  to  meters  per  second  in 
Tables  4.3  and  4.4.  Peak  particle  velocities  are  plotted  versus  acaled 
distance  in  Figures  4.9  and  4.10. 
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scaled  distance  (rtvVg-'  i/3) 


Figura  4.4  wava  valocity  as  a  function  of  acalad  diatanca  for  Tyndall 
baach  sand. 


Pleura  4. 
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100 
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SCALED  DISTANCE  (irVKQ  "  1/3) 


S  (cantinuad)  Wava  valocity  as  a  function  of  acalad  ( 
Ottawa  20-30  sand. 


6: 


PARTICLE  ACCELERATION  (9) 


Figur*  4.6  Typical  accel«cation-tijsc  history  and  intagratsd  valocity-tima 
history  curves. 
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PARTICLE  VELOCITT  Jg- 


T 


Tabla  4.3  ACCXLKRONKTSR  RESULTS  FOR  TYNDALL  BEACH  SAND 


Gage 

Peak 

Peak 

Charge 

Dist.  to 

Particle 

Particle 

Test 

Wt. 

Instr. 

Charge 

Accel. 

Velocity 

(gm) 

(cm) 

{m/s -2) 

(g-s)  (m/s) 

21.3-0-7 

1.031 

A1 

12.7 

121739 

0.22 

216 

(18.86) 

A2 

17.78 

45410 

0.12 

1.18 

A3 

20.32 

26280 

0.08 

0.78 

29.8-0-9 

0.35 

A1 

12.7 

52050 

0.11 

1,04 

(26.34) 

A2 

17.78 

A3 

20.32 

10824 

0.04 

0.35 

29.8-20-2 

0.35 

A1 

12.7 

68822 

0.12 

1.18 

(26.34) 

A2 

17.78 

24431 

0.06 

0.55 

A3 

20.32 

11215 

0.04 

0.34 

21.3-20-3 

1.031 

A1 

12.7 

159186 

0.24 

235 

(18.86) 

A2 

17.78 

56155 

0.12 

1.18 

A3 

20.32 

34814 

0.08 

0.78 

29.8-40-1 

0.35 

A1 

12.7 

48146 

0.1  i 

1.08 

(26.34) 

A2 

17.78 

9557 

0.04 

0.41 

.A3 

20.32 

8873 

0.03 

0  31 

21 .3-40-2 

1.031 

A1 

12.7 

164535 

0.25 

240 

(18.8o) 

A2 

17,78 

69447 

0.12 

1.20 

A3 

20,32 

46905 

0.08 

0.82 

29.8-60-1 

0.35 

A1 

12.7 

44821 

0.10 

0.94 

(26.34) 

A2 

17.78 

18379 

0.05 

0.45 

A3 

20.32 

16041 

0.04 

0.41 

21 .3-60-2 

1.031 

A1 

12.7 

141402 

0.27 

255 

(18.86) 

A2 

17.78 

59209 

0.12 

1.18 

A3 

20.32 

40128 

0.09 

0.36 

21.3-70-3 

1.031 

A1 

127 

152391 

0.23 

226 

(18.86) 

A2 

17.78 

49595 

0.11 

1.06 

A3 

20.32 

38957 

0.08 

0.82 

29.8-70-4 

0.3C 

A1 

12.7 

51327 

0.11 

1.06 

(26.34) 

A2 

17.78 

17813 

0.06 

0.54 

A3 

20.32 

13311 

0.04 

0.41 

PeU’entheses  indicate  g  level  at  instalment  depth. 
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Tabl.  4.4  ACCELEROKBTBR  RESULTS  FOR  OTTEWA  20-30  SAND 


Gan? 

Peak 

Peak 

Charge 

DIsl.  .0 

Particle 

Particle 

Test 

WL 

Instr. 

Charge 

Accel. 

Velocity 

(gm) 

(cm) 

(m/s ''2) 

(g-s)  (m/s) 

75.3-0-10 

0.35 

A3 

20.32 

2828.2 

0.02 

0.21 

(66.65) 

7S.3-r^  5 

0.35 

A1 

12.7 

35336.6 

0.10 

0.96 

(66.6i.^ 

A2 

17.78 

7502.7 

0.04 

0.34 

A3 

20.32 

3226.5 

0.03 

0.25 

75.3-40-5 

0.35 

A1 

1^7 

17638.4 

0.05 

0.47 

(66.65) 

A2 

17.78 

13798.7 

0.05 

0.50 

A3 

20.32 

8532.7 

0.04 

0.34 

75.3-60-S 

0.35 

A1 

1Z7 

28916.0 

0.10 

1.02 

(66.65) 

A2 

17.78 

9207.7 

0.04 

0.39 

A3 

20.32 

2781.1 

0.02 

0.22 

Parentheses  indicate  g  level  at  instrument  depth. 
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SCALED  DISTANCE  (m/kg  1/3) 


Figurs  4. 


8  Scaled  pealc  particle  acceleration  as  a  function  of  scaled 
distance  for  Ottawa  20-30  sand. 


s  ■  0\  H 


SCALED  QSTANCE  (m/kQ'‘l/3) 


Tifura  4.9  Paak  partiela  vaiocity  at  a  function  of  acalad  diatanea  for 
Tyndall  baaeh  aand. 


SECTION  V 


ANALYSIS  OF  RESULTS 

A.  PEAK  STRESS,  PEAK  PARTICLE  ACCELERATION,  AND  PEAK  PARTICLE  VELOCITY 

PREDICTION  EQUATIONS 

Prae-field  stresa  and  ground  motion  prediction  equations  from  Drake 
and  Little  (19R3)  ware  presented  in  Equations  (2.10)  through  (2.14)  of 
Section  II. A. 4.  Equations  were  developed  for  peak  stress,  peak  particle 
acceleration,  peak  particle  velocity,  peak  particle  displacement,  and 
impulse  and  are  expressed  as  an  equation  of  a  line  in  the  form  of  Equation 
(2,9) 

y  -  l)(x)*.  (2.9) 

If  the  constants  b  and  m  ara  known  for  a  data  sat,  a  response  equation  as 
a  function  of  scaled  distance  can  be  defined. 

Prediction  equations  for  peak  stress,  peak  particle  acceleration, 
and  peak  particle  velocity  on  Tyndall  beach  and  Ottawa  20-30  sands  are 
presented  in  this  section  in  the  form  of  Equation  (2.9).  Development  of 
response  equations  for  peak  particle  displacement  and  impulse  have  not 
been  pursued  and  will  not  be  analyzed  in  this  investigation. 

1.  Equations  for  Peak  Stress 

Plots  of  peak  stress  versus  scaled  distance  for  Tyndall  and 
Ottawa  sands  were  shown  in  Figures  4.2  and  4.3.  Values  of  b  and  m  are 
determined  from  these  plots  for  various  compaction  saturations  and  are 
listed  in  columns  one  (1)  through  three  (3)  of  Tables  5.1a  and  5.1b.  From 
Drake  and  Little's  (1983)  equation  for  stress  (Equation  (2.10)),  the 
constant  m  is  denoted  by  -n  and  the  constant  b  in  SI  units  is  eq'iivalent 


Table  5 -la  PEAK  STRESS,  PEAK  PARTICLE  ACCELERATION  AND  PEAK  PARTICLE  VELOCITY  EQUATION  CONSTANTS 
TYNDALL  BEACH  SAND 
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Not  determined  due  to  inadequate  data 


f0.049pCc(2.52)'*.  Paak  atraas  aquationa  davalopad  for  Tyndall  and 
Ottawa  aanda  ara  ahown  in  Tablaa  5.2  and  5.3. 

2.  Bquationa  for  Paak  Particle  Accalaration 

Prom  tha  acalad  paak  particle  acceleration  ragreaaiona  of 
Piguraa  4.7  and  4.8,  tha  conatanta  b  and  m  ara  obtained  for  both  Tyndall 
and  Ottawa  aanda  and  ara  liatad  in  columna  four  (4)  through  aix  (6)  of 
Tablaa  5.1a  and  5.1b.  Prom  Drake  and  Littla'a  (1983)  Equation  (2.12)  for 
peak  particle  accalaration,  m  ia  equivalent  to  -n-1  and  tha  constant  b  in 
SI  units  ia  fl26c«(2.52)*'.  Paak  particle  acceleration  aquationa  for 
Tyndall  and  Ottawa  sands  are  given  in  Tables  5.2  and  5.3. 

3.  Equations  for  Peak  Particle  Velocity 

From  the  peak  particle  velocity  ragreasions  of  Figures  4.9  and 
4.10,  the  constants  b  and  m  are  determined  for  both  Tyndall  and  Ottawa 
aanda  and  are  listed  in  columns  seven  (7)  through  nine  (9)  of  Tables  5. la 
and  5.1b.  The  constant  ro  is  equivalent  to  -n  and  b  in  SI  units  is 
£49(2.52)'^  from  Drake  and  Little's  (1983)  peak  particle  velocity  equation 
(Equation  (2.11)).  Particle  velocity  aquations  for  Tyndall  and  Ottawa 
sands  are  given  in  Tables  5.2  and  5.3. 

B.  ANALYSIS  OF  ATTENUATION  COEFFICIENTS 


The  elopee  of  tha  least  square  regression  lines  from  the  log-log 
data  for  peak  stress  (-n),  peak  particle  acceleration  (-n-1),  and  peak 
particle  velocity  (-n)  of  Table  5.1  share  the  common  variable  n  (Drake  and 
Little,  1983).  The  variable  n  was  defined  as  tha  attenuation  coefficient 
in  Section  II. A. 4.  Attenuation  coefficients  for  Tyndall  and  Ottawa  sands 
are  plotted  as  a  function  of  degree  of  saturation  during  compaction  in 
Figures  5.1(a)  and  (b) . 

Tyndall  sand  attenuation  coefficients  (Figure  5.1 (a;,  display  soma 
change  over  the  range  of  compaction  saturations  from  0  to  70  percent.  The 
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Tab(*  5.2a  STKSS  AW  G«OUW  NOTIOM  PRfDICTlOM  EOUATIOMS  OERIVB)  EMU  TtWALl  SAW  RECRESSIOH  DATA  (SI  UKITS). 
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26.41 


ATTENUATION  COEFFICIENT ,  n  ATTENUATION  COEFFICIENT , 


Figuc*  S.l  Atttnuation  eo«f f ieianta  «•  a  function  of  saturation  from 
ra^raaiiona  of  atraaa,  particia  accaiaration,  and  particla 
valoeity  and  from  Draka  and  Littla  (1983).  (a)  Tyndall  baac.^ 

aand.  (b)  Ottawa  20-30  aand. 
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attftnuation  coefflcianta  ara  in  fcha  aama  ranga  aa  tha  aaaumad  anvelopa  of 
attenuation  coefficianta  (rapraaantad  by  tha  dottad  line)  from  Dra)ce  and 
Xiittle'8  (1983)  denaa  and  loose  aand  data.  The  aaauned  envelope  waa 
conetructed  with  data  from  Table  2.1,  and  tha  dotted  line  waa  added  to 
repreeent  data  at  intermediate  eaturationa  which  waa  not  available  from 
Drake  and  Little* a  reaulta. 

Attenuation  coefficienta  for  Ottawa  aand  display  little  change  over 
the  range  of  compaction  aaturatione.  Theae  coefficienta  again  fall  in  the 
aaaumed  envelopea  of  attenuation  coefficienta  from  Drake  and  Little  (1983) 
( Figure  S. 1 (b) ) . 

C.  ANALYSIS  OF  INTERCEPTS  VALUES 

The  y-intercepte,  b,  of  the  least  aquare  regreeaione  for  peak 
atreas,  scaled  peak  particle  acceleration,  and  peak  particle  velocity  are 
given  in  Table  S.l.  intercepts  are  determined  at  a  scaled  distance  of  one 
for  both  SI  and  English  units  of  R/W*'^.  Intercept  values  in  SI  units  for 
Tyndall  and  Ottawa  sanda  are  plotted  in  Figures  5.2(a),  (b),  and  (c)  as 
a  function  of  compaction  saturation.  Intercept  values  for  peak  stress  and 
scaled  peak  particle  acceleration  are  lowest  at  0  and  53  percent 
saturations  and  are  at  a  maximum  at  35  percent  saturation.  Intercepts  for 
peak  particle  velocity  are  lowest  at  0  and  70  percent  saturations  and 
again  are  at  a  maximum  at  35  percent  saturation.  Minimum  and  maximum 
intercept  values  for  both  peak  stress  and  scaled  peak  particle 
acceleration  differ  by  a  factor  of  two  in  some  instances. 

Peak  stress  intercepts  for  Ottawa  sand  follow  a  similar  trend  to 
stress  intercepts  of  Tyndall  sand.  Due  to  a  lack  of  accelerometer  data 
for  Ottawa  sand,  representative  trends  could  not  be  determined  for  scaled 
peak  particle  acceleration  and  peak  particle  velocity  intercepts.  This 
data  is  represented  by  a  dotted  line. 

Intercepts  for  Tyndall  and  Ottawa  sands  in  English  units  are  shown 
in  Figures  5.3(a),  (b),  and  (c)  as  well  and  exhibit  the  same  trends  as 
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P£=CSNT  SATURATTCN 


Tigura  S.2  Intarcapt  valua*,  b,  for  Tyndall  baach  and  Ottawa  20>30  sands 
in  SX  units.  (a)  Paak  strass.  (b)  Scalad  paak  particla 
aecalaration.  (e)  Paak  particla  valocity. 
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Mgur*  5. 


3  Intareapt  valusa,  b,  fer  Tyndall  batch  and  Ottawa  20-30  sanda 
in  Xngliah  ur.ita.  (a)  Peak  atraaa.  (b)  Scalad  paak  partiela 
aeealaration.  <e)  paak  partiela  velocity. 
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those  in  SI  units.  English  unit  intercept  values  can  be  obtained  from  SI 
unit  regression  lines  at  a  scaled  distance  of  0.4  m/kg''^. 

Wave  velocity,  stress  transmission  ratio,  peak  stress,  and  shear 
modulus  for  sands  compacted  moist  (Ross  et  al.,  1988;  Pierce  et  al.,  1989; 
Charlie  et  al.,  1990;  and  Wu  et  al.,  1984)  display  a  remarkably  similar 
trend  to  the  intercept  values  of  Figures  5.2  and  5.3  where  maximum  values 
occur  at  intermediate  saturations  (Section  Il.B.l.). 

Bishop  and  Blight  (1963)  suggested  that  an  increase  in  effective 
stress  in  unsaturated  soils  is  due  to  an  increase  in  the  soils  capillary 
pressure.  The  effective  stress  increase  is  the  product  of  the  soil 
saturation,  S,  and  matric  suction,  (u,-u„) .  Charlie  et  al.  (1990) 
observed  the  effects  of  matric  suction  on  capillary  stresses  for  silica 
50/80  sand.  The  result  was  a  capillary  stress  of  0  at  a  saturation  of  0 
and  100  percent  and  a  maximum  capillary  stress  at  approximately  38  percent 
saturation  (Figure  5.4).  Charlie  et  al.  (1990)  conducted  one  dimensional 
quasi-static  compression  tests  at  a  constant  strain  rate  on  dense  silica 
sand  to  determine  the  constrained  modulus,  M,  as  a  function  of  saturation. 
Modulus  values  for  the  stress-strain  relationships  increased  up  to  34 
percent  saturation  and  then  decreased  with  increasing  saturation  (Figure 
5.5).  Modulus  values  followed  the  same  trend  as  Figure  5.4  for  the 
product  of  matric  suction  and  saturation,  and  thus  demonstrated  a 
relationship  between  the  degree  of  saturation  during  compaction  and  soil 
modulus. 

Comparing  the  results  of  Figure  5.5  for  Silica  50/80  sand  and  the 
intercept  results  of  Figures  5.2  and  5.3,  we  conclude  that  stress  wave 
propagation  is  a  function  of  modulus  as  well  as  compactive  saturation. 

D.  ANALYSIS  OF  SEISMIC  WAVE  VELOCITY  AND  CONSTRAINED  MODULUS 

At  large  distances  from  a  blast  soil  strains  are  small  and  are  part 
of  the  soils  elastic  range  (strains  less  than  10'*%).  Small  strain  wave 
velocity,  or  seismic  velocity,  is  a  basic  index  of  the  dynamic  response 
of  laoils. 
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Figure  S.4  Increase  in  effective  streae  for  50/80  Silica  sand  due  to 
saturation  and  matrix  suction  (From  Charlie  et  al.,  1990). 
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Figure  5.5  Quasi-static  stress  strain  curves  for  50/80  Silica  sand  (From 
Charlie  at  al..  1990). 


A  plot  of  wave  velocity  versus  scaled  distance  for  dry  Tyndall  beach 
sand  is  shown  in  Figure  5.6(a).  The  seismic  velocity  is  approached  as  the 
scaled  distance  increases  and  the  wave  velocity  becomes  relatively 
constant.  An  approximate  range  of  seismic  velocity  can  be  determined  from 
a  log-log  plot  of  the  data  at  a  scaled  distance  of  2.9  m/kg''’  (Figure 
5.6(b)).  The  elastic  range  for  these  sands  is  calculated  to  be  at  scaled 
distances  greater  than  10  m/kg'^. 

Seismic  wave  velocity  approximations  for  Tyndall  sand  are  shown  in 
column  one  (1)  of  Table  5. 4a  and  for  Ottawa  sand  are  shown  in  column  one 
(1)  of  Table  5.4b.  Seismic  velocity  as  a  function  of  saturation  is 
plotted  in  Figures  5.7(a)  and  (b)  . 

Seismic  velocity  taken  from  Table  2.1  by  Drake  and  Little  (1983)  for 
dense  and  loose  sands  are  shown  in  column  two  (2)  of  Tables  5.4a  and  5.4b. 
Seismic  velocities  at  intermittent  saturations  are  not  reported  by  Drake 
and  Little  (1983)  and  are  represented  by  an  envelope  of  assumed  data 
(Figures  5.7(a)  and  (b)). 

Peak  stress  and  peak  particle  velocity  are  related  to  one  another 

through  Equation  (2.3).  Seismic  velocity  can  be  expressed  in  terms  of 

peak  stress  and  peak  particle  velocity  by  rearranging  the  terms  of 

equation  (2.3)  Using  this  relation,  seismic  velocities  for  Tyndall  sand 

(5.1) 


are  determined  from  peak  stresses  and  peak  particle  velocities  of  Figures 
4.2  and  4.9  and  are  shown  in  column  three  (3)  of  Table  5.4a.  Peak 
stresses  and  peak  particle  velocities  for  Ottawa  sand  are  determined  from 
Figures  4.3  and  4.10,  and  computed  seismic  velocities  are  given  in  column 
three  (3)  of  Table  5.4b.  As  shown  in  Figures  5.7(a)  and  (b),  seismic 
velocities  computed  using  Equation  (5.1)  are  comparable  to  the  range  of 
measured  wave  velocities  and  seismic  velocities  determined  by  Drake  and 
Little  (1983). 

Peak  stress  and  peak  particle  velocities  of  Equation  (5.1)  are 
assumed  to  be  in  the  elastic  range  of  the  soil,  therefore  resulting  in 
seismic  velocities  of  the  elastic  range  as  well.  Peak  stress  and  peak 
particle  velocity  values  for  Tyndall  and  Ottawa  sands  in  this 
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SCALED  DISTANCE  (m/kg''  1/3) 
(b) 


Figure  5.6  Wave  velocity  versus  scaled  distance  for  dry  Tyndall  beach 
sand.  (a)  Normal  scale.  (b)  Log-log  scale  approximate 
siesmic  velocities  at  scaled  distance  of  2.9  m/kg'^. 


Tablt  S.4s  COMPARISON  OF  WAVE  VELOCITIES  FOR  TYNDALL  BEACH  SAND. 

Wave  Valocity  <nv/s) 


N«atur«d  Drake  A  Littia 

SATURATION  Wave  Velocity  (1983)  Eq.  S.1  Eq.  S.8 

<1)  (2)  (3>  (4) 


0 

200-260 

207 

384 

17 

2S0-360 

314 

20 

275-396* 

31 

155* 

35 

210-290 

369 

53 

220-370 

305 

70 

240-400 

220 

74 

300* 

76 

152-183* 

*  Determined  from  dense  poorly  graded  sand. 

*  Determined  from  loose  poorly  graded  sand. 


Table  5.4b  COMPARISON  OF  WAVE  VELOCITIES  FOR  OTTAWA  SAND. 


Wave  Velocity  (n/s) 

Measured  Drake  (  Little 

SATURATION  Wave  Velocity  (1983)  Eq.  5.1  Eq.  5.8 

(1)  (2)  (3)  (4) 

ISOiTiO  S  5W~ 

220  275- 396*  174 

183* 

230-260  94 

200-240  146 

300* 

152-183* 


*  Not  determined  due  to  inadequate  data. 

*  Determined  from  dense  poorly  graded  sand. 

*  Determined  from  loose  poorly  graded  sand. 


20 

31 

40 

60 

74 

76 
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WAVE  VELOCITY  (m/s)  WAVE  VELOCITY  (m/s) 


Figure  5.7  Meaaured  wave  velocitiee  from  centrifuge  data  at  a  scaled 
range  of  2.9  m/kg''’«  and  eeiemic  velocitiee  from  Drake  and 
Litrle  (1983),  and  comt>u’<.ed  from  Equationa  (S.l)  and  (5.8). 
(a)  Tyndall  beach  sand.  (b)  Ottawa  20-30  sand. 


inveetlgation  are  not  available  from  the  elastic  range.  However,  values 
can  be  obtained  from  a  scaled  distance  of  2.9  m/kg''*  and  will  be  used 
herein.  The  density,  p,  for  use  in  Equation  (5.1)  for  moist  specimens  is 
assumed  constant  as  a  result  of  desaturation  to  the  residual  saturation 
(APPENDIX  C). 

Seismic  velocity  and  constrained  modulus  are  related  to  one  another 
through  Equation  (2.4).  Rearrangetnent  of  Equation  (2.4)  allows  one  to 
solve  for  M  in  terms  of  c, 

M  ■  pc|.  (5.2) 

Again,  the  mass  density  of  moist  centrifuge  specimens  is  assumed  constant 
due  to  desaturation.  Modulus,  therefore,  is  directly  proportional  to  the 
square  of  the  seismic  velocity.  Values  of  constrained  modulus  for  Tyndall 
sand  are  computed  from  seismic  velocities  given  in  column  one  (1)  of  Table 
5.4a,  and  are  listed  in  column  one  (1)  of  Table  S.Sa.  The  proportional 
relationship  of  seismic  velocity  and  modulus  is  evident  as  they  both 
display  the  same  unique  trend.  Modulus  values  are  also  determined  from 
measured  seismic  velocities  (column  one  (1),  Table  5.4a)  and  listed  in 
column  two  (2)  of  Table  5. 5a. 

Values  of  constrained  modulus  for  Tyndall  sand  are  plotted  versus 
saturation  in  Figure  5. 8(a).  For  comparison,  assumed  modulus  envelopes 
for  dense  and  loose  sands  by  Dra)(«>  and  Little  (1983)  are  included  in 
Figure  5.8(a).  Values  of  constrained  modulus  for  Ottawa  sand  are  computed 
similarly  to  those  of  Tyndall  sand  and  are  presented  along  with  values 
from  Drake  and  Little  (1983)  in  Table  5.5b.  A  plot  of  modulus  versus 
degree  of  saturation  for  Ottawa  sand  is  shown  in  Figure  5.8(b). 

One  dimensional  quasi-static  compression  teats  at  constant  strain 
rate  wore  conducted  on  Tyndall  sand  to  experimentally  determine  the 
constrained  modulus  at  large  strains.  These  results  are  given  in  column 
five  (5)  of  Tadile  5.Sa  along  with  the  large  strain  modulus  results  of 
Charlie  et  al.  (1990)  from  Section  V.B.  (column  (4),  Table  5.5a). 
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Tablt  S.4a  COHPAKISON  OF  UAVE  VELOCITIES  FOA  TYNDALL  BEACH  SAND. 


Uavc  Valocity  (m/s) 


SATURATION 

Measured 
Wave  Velocity 
(1) 

Drake  (  Little 
(1983) 

(2) 

Eq.  5.1 
(3) 

Eq.  S.a 
(4) 

0 

200-260 

207 

384 

17 

250-360 

314 

20 

275-396* 

31 

183“ 

35 

210-290 

369 

S3 

220-370 

305 

70 

240-400 

220 

74 

300* 

76 

152-183“ 

*  Dctsminad  from  da.tse  poorly  graded  sand. 
'  Determined  from  loose  poorly  graded  sand. 


Table  5. 4b  CamPARISON  OF  WAVE  VELOCITIES  FOA  OTTAWA  SAND. 

Wave  Velocity  <m/s) 


Measured  Drake  (  Little 

SATURATION  Wave  Velocity  (1983)  Eq.  5.1  Eq.  5.8 

(1)  (2)  (3)  (4) 


0 

130-210 

• 

525 

20 

220 

275-396* 

174 

31 

183* 

40 

230-260 

94 

60 

200-240 

146 

74 

300* 

76 

152-183“ 

*  Not  detersiined  due  to  inadequate  data. 

*  Determirwd  from  derwe  poorly  graded  sand. 

*  Determined  from  loose  poorly  graded  sand. 
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MODULUS  (MPa)  MODULUS  (MPa) 


FEnCENT  SATURATION 

(b) 


FiTura  5.8  Constrainad  modulus  as  a  function  of  saturation  cemputad  from 
Equation  15 .2  I  and  takan  from  Draka  and  Littla  (1983  1 
Tyndall  b'aach  sand,  (b)  Ottawa  20-30  sand. 
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E. 


ANALYSIS  OF  ACOUSTIC  IMPEDANCE 


The  acoustic  impsdsncSf  pc,,  is  a  isaasure  of  a  soils  ability  to 
transmit  stress  waves.  Using  Drake  and  Little's  (1983)  peak  stress 
equation  (Equation  (2.10)),  one  can  back  calculate  pc,  as  a  function  of 
stress  snd  scaled  distance 


PCc 


P 


f  0.049  [: 

i.52— ^1 

1^1/1  j 

(5.3) 


Acoustic  impedance  at  a  scaled  distance  of  2.9  m/kg*'’  is  given  in  column 
one  (1)  of  Tables  S.6a  and  S.6b  for  Tyndall  and  Ottawa  sands.  Back 
calculating  pc,  using  the  above  equation  assumes  that  the  constant 
fO.  049pc,(  2 . 52  )'*,  empirically  determined  by  Drake  and  Little  (1983),  is 
correct . 

Acoustic  impedance  can  also  be  computed  as  a  function  of  peak  stress 
and  peak  particle  velocity  in  a  form  similar  to  Equation  (5.1). 


pCc  = 


P 

V' 


(5.4) 


Values  of  pc,  from  Equation  (5.4)  are  given  in  column  two  (2)  of  Tables 
5. 6a  and  5.6b  for  Tyndall  and  Ottawa  sands.  Plots  of  pc,  for  Tyndall  and 
Ottawa  sands  are  shown  in  Figures  5.9(a)  and  (b)  and  include  assumed 
envelopes  of  pc,  observed  by  Drake  and  Little  (1983). 


F.  ANALYSIS  OF  PEAK  STRESS 


Envelopes  of  peak  stress  data  taken  from  Drake  and  Little  (1983)  for 
very  loose  dry  sand,  loam  loess  and  dry  sand,  and  very  dense  sand  are 
shown  in  Figure  5.10  along  with  peak  stress  data  for  dry  dense  and  medium 
dense  Tyndall  and  Ottawa  sands.  Both  attenuation  and  magnitude  of  peak 
stresses  for  the  two  sets  of  data  are  similar.  Centrifuge  data  also 
extends  over  a  scaled  distance  similar  tc  that  of  Drake  ana  Little  (also 
see  Figures  4.2  and  4.3). 
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Table  5.6a  COMPARISON  OF  ACOUSTIC  IMPEDANCE  FOR  TYNDALL  BEACH  SAND 


ACOUSTIC  IMPEDANCE  X  10’  (kg/m’-s) 

Drake  &  Little 


SATURATION 

Eg.  5.3 
<1) 

Eg.  5.4 

(2) 

(1983) 

(3) 

0 

272 

315 

17 

227 

492 

20 

566* 

31 

262'’ 

35 

232 

579 

S3 

225 

478 

70 

285 

345 

74 

498' 

76 

283-339*’ 

‘  Determined  from  dense  poorly  graded  sand. 
^  Determined  from  loose  poorly  graded  sand. 


Table  S.6b  COMPARISON  OF  ACOUSTIC  IMPEDANCE  FOR  OTTAWA  SAND. 


ACOUSTIC  IMPEDANCE  X  lO’  (kg/m’-s) 


SATURATION 

Eg.  5.3 
(1) 

Eg.  5.4 
(2) 

Drake  6  Little 
(1983) 

(3) 

0 

290 

• 

20 

447 

283 

566* 

31 

262*’ 

40 

278 

154 

60 

226 

238 

74 

498‘ 

76 

283-339*’ 

*  Not  determined 

due  to  inadeguate  data. 

*  Determined 

from 

danse  poorly  graded 

sand. 

*  Determined 

from 

loose  poorly  graded 

sand. 
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ACOUSTIC  IMPEDANCE  (kg/m  2-s)  ACOUSTIC  IMPEDANCE  (kg/m  2-s) 


■  Equation  S.3 
■I-  Equation  5.4 
^  Drake  and  Little  (1S83) 


STRESS  (kPa)  PEAK  STRESS  (hPa) 


Plots  of  p^ak  atrssa  as  a  function  of  peak  particle  velocity  for 


Tyndall  beach  sand  are  shown  in  APPENDIX  D. 


G.  ANALYSIS  OF  UNSATURATED  SOIL  MECHANICS  THEORY 

Unsaturated  soil  mechanics  theory  is  based  on  the  linear  behavior 
of  soil  properties,  such  as  shear  strength  and  strain,  between  dry  and 
saturated  states.  For  instance,  unsaturated  soil  mechanics  theory  by 
Fredlund  (1985),  described  in  Section  II. B. 3,  assumes  that  shear  strength 
of  soil  varies  linearly  with  (o-u.)  and  (u,-u.)  (Figure  2.10).  Similarly, 
Fredlund's  (1985)  equation  for  soil  strain 

£  ■  ffii'  d(o-u^)  cHu^-uJ  (2.22) 

utilizes  the  same  stress  components  as  does  shear  ((o-u.)  and  (u,**u^)),  and 
given  that  the  compressibility  factors  m,*  and  it;*  are  constant,  then 
resulting  soil  strains  vary  linearly  as  well. 

Pierce  (1989)  derived  an  equation  to  predict  modulus  for  dense  Eglin 
sand  by  utilizing  the  experimental  results  from  split~Hopkinson  pressure 
bar  testing  at  various  values  of  (o-u,)  and  (u,-u«).  Pierce  based  Equation 
(5.S)  on  Fredlund's  (1985)  Equation  (2.22)  for  strain. 

M-  (1.510  -  0.016(u,-uJ  ♦  0.0004  (o-u.)  )  X"  10*  (5.5) 

()iPa).  Pierce  observed  that  the  modulus  of  Eglin  sand  was  relatively 
unaffected  by  changes  in  pore  water  pressures  but  was  affected  by  changes 
in  total  strees.  Tests  specimens  were  compacted  dry,  saturated,  and  then 
the  water  content  was  adjusted  using  the  pressure  plate  method.  Pierce's 
results  are  supported  by  Fredlund's  theory  that  pore  water  pressure  has 
little  influence  on  changing  the  modulus  of  the  sand. 

Results  reported  by  Charlie  et  al.  (1990)  for  sands  compacted  moist 
are  not  predicted  by  Fredlund's  (1985)  theory.  As  demonstrated  in  Figure 
5.5,  modulus  is  related  to  the  saturation  at  the  time  of  compaction  and 
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displays  non-linear  trends  over  the  range  of  saturations  from  0  to  66 
percent.  Since  modulus  is  directly  related  to  soil  strength,  soil 
strength  is  assumed  to  vary  non-linearly  with  the  degree  of  saturation  at 
the  time  of  compaction  as  well.  This  would  infer  that  soil  properties, 
such  as  shear  strength  and  soil  strain,  also  vary  non-linearly  with 
modulus  and  compactive  saturation. 

Studies  by  Hardin  and  Black  (1968)  on  various  quantities  which 
influence  shear  modulus,  G^,  of  sand  have  shown  that  for  shearing  strain 
amplitudes  lass  than  10**%,  Q,  was  essentially  independent  of  each  of  tne 
variables  except  effective  stress,  o',  and  void  ratio,  a.  Saturation  had 
little  effect  on  G,.  Hardin  and  Richart  (1963)  developed  an  empirical 
expression  to  relate  G.,  e,  and  o'  for  round-grained  soils 

^  2l£(2^_17^^^  (kPa)  .  (5.6) 


By  applying  concepts  of  the  theory  of  elasticity  and  Equations  (2.6)  and 
(2.7),  shear  modulus  can  be  related  to  constrained  modulus,  M,  in  the  form 


M  « 


2(l-M)go 

l-2p 


(5.7) 


Substituting  constrained  modulus  into  Equation  (2.4)  yields  seismic 
velocities  as  a  function  of  e  and  o'  in  the  form 


p  (l-2n)  (1+e) 


(5.8) 


Although  Equation  (5.8)  is  easily  derived,  its  validity  for  moist  soils 
is  questionable.  The  constrained  modulus  of  Equation  (5.7)  is  defined  as 
a  function  of  shsar  modulus,  and  shear  modulus  for  sand  is  experimentally 
determined  dry.  The  use  of  Equation  (5.8)  to  derive  seismic  velocities 
becomes  further  complicated  when  determining  effective  stress  of  an 
unsaturatsd  soil.  One  must  then  rely  on  unsaturated  theory  such  as 
Equation  (2.17)  from  Bishop  (1960). 

Seismic  velocities  for  dry  Tyndall  and  Ottawa  sands  using  Equation 
(5.8)  are  given  in  column  four  (4)  of  Tables  5.4a  and  5.4b. 
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SECTION  VI 


SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 


A>  SUMMARY 


Th«  survivability  of  a  buried  structure  froa  a  nearby  blast  loading 
is  dependent  on  the  dynamic  soil  properties  of  the  backfill  material 
around  and  over  the  structure.  Moisture  in  the  backfill  material  at  the 
time  of  compaction  has  been  shown  to  affect  the  soil  stiffness,  therefore, 
influencing  the  soils  stress  transmission  ability.  This  report  presents 
the  results  of  a  series  of  systematic  explosive  tests,  in  representative 
cohesionless  backfill  material,  conducted  to  determine  the  effect  of 
saturation  at  the  time  of  compaction  on  blast  induced  stress  wave 
propagation  and  ground  motion  response.  Tyndall  beach  and  Ottawa  20/30 
sands  were  utilized  for  the  backfill  material. 

Explosive  testing  was  conducted  on  the  centrifuge  in  soil  models  of 
1/18.9  and  1/26.3  scale  and  were  accelerated  to  16.9  and  26.3  g's 
respectively.  With  the  use  of  scaling  laws,  detonators  having  1031  mg  and 
350  mg  of  PBX  9407  were  used  to  simulate  6.9  kg  and  6.4  kg  (7.8  kg  and  7.3 
kg  THT  equivalent)  prototype  charges  respectively.  Charges  were  buried 
to  simulate  contained  bombs  at  a  depth  of  burial  of  1.4  meters  resulting 
in  a  ground  shock  coupling  factor  of  0.71  m/kg''^  (1.8  ft/lb'^).  By  using 
models  to  simulate  prototype  blast  events,  tasting  could  be  conducted  at 
a  fraction  of  the  money  and  time  required  to  perform  large  scale  blast 
testing. 

Model  and  prototype  performance  are  related  using  a  set  of 
parameters  known  as  scaling  laws.  The  Buckingham  Pi  Theory  was  used  to 
develop  these  scaling  laws  through  the  principles  of  dimensional  analysis. 
By  using  two  model  charge  masses  to  represent  a  single  prototype  TNT 
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charge  maea  of  approximately  7.5  kg,  the  principle  of  "modeling  of  modele" 
was  implemented  and  provided  a  check  on  the  consistency  of  the  scaling. 

Models  were  compacted  moist  to  target  saturation  levels  and  a 
constant  void  ratio  using  a  method  of  vibration.  Compaction  by  vibration 
was  chosen  to  be  the  most  efficient  means  for  preparing  soil  models 
because  of  the  ability  to  reproduce  uniform  specimens  and  oppose  the 
forces  of  capillarity. 

From  blast  response  data  for  Tyndall  and  Ottawa  sands,  empirical 
equations  for  peak  stress,  peak  particle  acceleration,  and  peak  particle 
velocity  were  developed  as  a  function  of  compaction  saturation. 
Attenuation  coef f icients,  wave  velocities,  constrained  modulus,  and 
acoustic  impedance  for  Tyndall  and  Ottawa  sands  were  compared  with  data 
from  Drake  and  Little  (1983).  Intercept  values  at  scaled  dictances  of 
1  m/kg''’  and  1  ft/lb''^  were  analyzed  for  peak  stress,  scaled  peak  particle 
acceleration,  and  peak  particle  velocity. 

6.  CONCLUSIONS 


1.  Stress  and  Oround  Motion  Empirical  Equations 

Stress  and  ground  motion  equations  are  obtained  from  least 
square  regression  lines  of  soil  response  data  as  a  function  of  scaled 
distance.  These  equations  are  defined  by  two  components:  the  slope  of  the 
regression  line  (ro)  or  attenuation,  and  the  y-intercept  of  the  line  (b) 
at  a  scaled  distance  of  one.  Slope  and  intercept  values  as  a  function  of 
compaction  saturation  for  Tyndall  and  Ottawa  sands  are  listed  in  Tables 
5.1a  and  5.1b.  Prediction  equations  developed  from  these  constants  for 
peak  stress,  peak  particle  acceleration  and  peak  particle  velocity  ere 
shown  in  Tables  5.2  and  5.3. 

Peak  stress,  scaled  peak  particle  acceleration,  and  peak  particle 
velocity  intercept  values  for  both  Tyndall  and  Ottawa  sands  are  strongly 
influenced  by  compaction  saturation.  Tyndall  sand  intercepts  for  peak 
stress  and  scaled  peak  particle  acceleration  are  lowest  at  0  and  53 
percent  saturations  and  are  at  a  maximum  at  35  percent  saturation  (Figures 


102 


S.2(a)  and  (b)).  Intercepta  for  paak  particla  valocity  are  loweat  at  0 
and  70  percent  aaturationa  and  again  are  at  a  maximum  at  35  percent 
saturation  (Figure  5.‘>(c)).  Peak  atress  intercept  values  for  Ottawa  sand 
follow  a  similar  trend  to  stress  intercept  values  of  Tyndall  sand  (Figure 
S.3(c)).  Minimum  and  maximum  intercept  values  for  Tyndall  scaled  peak 
particle  accelerations  and  Ottawa  peak  stresses  differ  by  a  factor  of  two. 
Because  of  Limited  data  obtained  for  Ottawa  sand,  intercept  trends  could 
not  be  observed  for  scaled  peak  particle  acceleration  and  peak  particle 
velocity. 


2.  Comparison  with  Split-Hopkinson  Pressure  Bar  Tests 

Wave  velocity,  constrained  modulus,  and  acoustic  impedance 
were  computed  from  response  data  for  Tyndall  sand.  Results  display 
minimum  values  at  0  and  70  percent  saturations  and  are  maximum  at  the  mid 
saturations  (Figures  5.7,  5.8  and  5.9).  Similar  trends  were  observed  for 
Tyndall  and  Ottawa  sand  intercept  values.  This  trending  pattern  closely 
compares  to  the  trending  pattern  for  wave  velocity  and  stress  transmission 
ratio  observed  by  Ross  et  al.  (1988),  Pierce  et  al.  (1989),  and  Charlie 
et  al.  (1990)  on  the  SHPB  (Section  II.B.l.).  Centrifuge  results  thus 
indicate  that  an  increase  in  blast  induced  stress  transmission  occurs  for 
soils  compacted  over  the  mid  saturation  range. 

3.  C-  .-i  •  ,ith  Drake  and  Little  (j.933)  Blast  Tasting 

Both  attenuation  and  magnitude  of  peak  stress  as  a  function 
of  scaled  distance  for  dry  Tyndall  and  Ottawa  sands  compare  closely  to 
peak  stress  data  from  Drake  and  Little  (1983)  (Figure  5.10).  Attenuation 
coefficients,  wave  velocity,  constrained  modulus,  and  acoustic  impedance 
determined  from  centrifuge  modeling  closely  approximate  proximity  to  those 
obtained  by  Drake  and  Little  (Figures  5.1,  5.7,  5.8,  and  5.9).  This 
demonstrates  the  centrifuges  ability  to  model  large  scale  stressba  as  a 
result  of  blast  loading.  Verification  of  prototype  explosive  events  is 
currently  being  conducted  by  Villano  (1992)  and  Dowden  (1992). 
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4.  Comparison  with  Unsaturatsd  Soil  Mechanics  Theory 

Finally,  results  by  Ross  et  al.  (1988),  Pierce  et  al.  (1989), 
and  Charlie  et  al.  (1990),  along  with  the  results  presented  in  this 
investigation,  do  not  agree  with  unsaturated  soil  mechanic  theory  by 
Fredlund  (1985).  As  demonstrated  in  this  investigation,  the  non-linear 
behavior  of  constrained  modulus  as  a  function  of  compactive  saturation 
does  not  follow  the  linear  shear  stress  and  matric  suction  relationship 
proposed  by  Fredlund  (Figure  2.10).  Centrifuge  results  indicate  that 
greater  shear  stresses  are  developed  in  soils  compacted  in  the  mid 
saturations . 

C.  RECOMMENDATIONS  FOR  DESIGN 


Explosive  test  results  for  Tyndall  beach  and  Ottawa  20-30  sands 
compacted  moist  (saturations  ranging  from  approximately  20  to  50  percent) 
indicate  that  Equations  (2.10)  to  (2.12),  developed  by  Drake  and  Little 
(1983),  underpredict  peak  stress,  peak  particle  velocity  and  peak  particle 
acceleration  by  a  factor  of  about  two.  Therefore,  for  cohesionless 
backfill  sands  compacted  moist,  it  is  recommended  that  the  predicted 
values  for  peak  stress,  peak  particle  velocity  and  peak  particle 
acceleration  be  doubled  for  design  purposes. 
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APPENDIX  A 


DERIVATION  OP  SCALING  LAWS 


The  Buckingham  Pi  theory  acta  ab  a  meane  to  relate  functions 
expressed  in  terms  of  dimensional  units  and  functions  of  non-dimensional 
units.  A  non-dimensional  quantity  is  defined  as  any  quantity  formed  in 
such  a  way  that  all  the  units  identically  cancel.  In  evaluating  models, 
it  is  necessary  to  define  a  function  relating  all  the  parameters  of  the 
modeled  system.  Buckingham  (1915)  proposed  a  theory  which  eliminates  the 
need  to  know  the  function,  allowing  direct  similitude  through 
non-dimensional  m  terms. 

In  any  given  physical  system  there  exists  one  or  more  dependent 
parameters  which  in  turn  are  a  function  of  independent  parameters. 
Independent  parameters  are  those  quantities  which  are  necessary  to  fix 
location  inside  a  given  problem,  whereas  dependent  parameters  are 
dependent  upon  the  location.  If  qi  is  a  dependent  parameter  and  q^,  q,, 
...  q,g  are  independent,  we  may  write 

<7i  ■  f'(<7a.<7a.  •  •  •  ,  G,)  •  (A.l) 


Equation  (A.l)  can  be  mathematically  expressed  as 

f<Gj.<5rarG3.  •  •  .  .<7,)  "0  (A. 2) 

and  parameters  can  be  combined  to  form  non-dimensional  tr  terms  such  that 

f (Rj , X],  .  .  . , Ti„)  ■  0  (A.3) 

The  number  of  n  terms,  n,  is  related  by 
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n  ^  m  ~  k 


{K.A) 


wher*  m  is  the  number  of  q'e,  end  k,  from  Buckingham'e  (1915)  original 
formulation  of  the  theorem,  ia  equal  to  the  minimum  number  of  independent 
dimeneione  required  to  conatruct  the  dimensions  of  all  other  parameters 
qii  qi,  •••  q..  In  conclusion,  if  all  of  the  relevant  parameters  are 
included  in  Equation  (A. 2),  then  similitude  between  model  and  prototype 
can  be  achieved  with  the  dimenaionlesa  n  terns  such  that 

“  *prototyp«*  (A*  5) 

Derivation  of  the  scaling  laws  of  Table  2.2  are  based  on  Equation 

(A. 5).  Ab  an  example,  using  the  dimensionless  tr  term,  Lgp/E,  we  can 

assume  where  m  represents  the  model,  p  represents  the  prototype  and 

(A. 6) 


L  ■  linear  dimension; 
g  «  gravitational  force; 
p  «  mass  density;  and 
E  «  modulus  of  elasticity. 

The  majority  of  scaling  literature  assumes  that  materials  used  in  the 
prototype  are  the  same  as  what  is  used  to  construct  the  model,  or  at  least 
that  the  material  properties  will  be  kept  constant  (Bradley  at  al.,  1984). 
Therefore,  E.  ■  and  ■  p^,  and  Equation  (A. 6)  becomes 

-  I-pSTp.  (A. 7) 

Since  the  prototype  is  assumed  to  be  at  1  g,  g,  •  1,  and  setting  g.  ■  N, 
we  arrive  at  the  first  scaling  law 

Lp  -  LJf.  (A. 8) 

The  length  term,  L,  and  the  gravity  term,  N,  are  used  to  derive  the  other 
scaling  laws  of  Table  2.2.  A  complete  derivation  of  Table  2.2  is 
conducted  by  Bradley  et  al.  (1984). 
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APPENDIX  B 


INSTRUMENTATION  CALIBRATION 

A.  STRESS  GAGE  CALIBRATION  RESULTS 

Quaei-atatic  atrass-voltage  ralationshipa  for  the  carbon  atreaa 
gages  were  experimentally  determined.  This  procedure  consisted  of 
vertically  loading  eight  4-inch  columns  of  Tyndall  Beach  sand  in 
increments  of  2.2  kN  (500  lb).  The  sand  was  confined  in  rigid  stainless 
steel  cylinders  (Figure  B.l)  with  a  carbon  gaga  placed  at  a  depth  of  2.54 
cm,  A  calibrated  load  cell  attached  to  the  loading  machine  was  used  for 
direct  load  readings  while  corresponding  carbon  resistance  change  was 
measured  with  a  digital  volt  meter  after  being  amplified  by  n  Bctron 
signal  conditioner.  Voltages  for  amplification  gain  were  corrected. 

For  pressures  up  to  19  MPa,  linear  stress-voltage  relationships 
resulted  (Figure  B.2).  Eight  resistors  were  calibrated  in  this  manor  and 
a  mean  regression  slope  of  190.7  NPa/volt  was  used  for  the  calibration  of 
all  gags  resistors  used  in  this  research  effort  (Figure  B.3}. 

B.  ACCELEROMETER  CALIBRATION  RESULTS 

Calibration  of  accelerometers  for  this  research  was  conducted  by  the 
manufacturer  (ENDEVCO).  Calibration  .urves  for  the  accelerometers  and 
their  sensitivities  are  shown  in  Figure  B.4.  All  calibrations  are 
traceable  to  the  National  Bureau  of  Standards. 


Ill 


pr-1. 


e  B.2  (continued)  Stress-voltage  relationships  and  regressions  for 

carbon  stress  gages. 


STRESS  (MPa) 


n^urtt  b.3  Si:r«ss  gag*  calibration  •nvalopc  and  maan  ragraaaion. 
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SHOCK  CALIBRATION 


TIME  J0O_uS«c/Dlv  SHOCK  AMPUTUOEi2iS_9 '  »PK 


SHOCK  calibration 


TIME  leq  uSae/Piv  SHOCK  AMPLITUDE  o'lQk 


SHOCK  CALIBRATION 


TIME  laa  ua«c/Dlv  shock  amplitude  *3Z3  g'lOK 


Figurt  B.4  Calibration  curves  Cor  Accelerometers 
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APPENDIX  C 


XFPBCTS  OP  DBSATURATION  ON  SPECIMEN  DENSITY 

Just  as  soil  dspth  scales  with  acceleration  of  the  centrifuge 
(Section  III.6.3.),  height  of  capillary  rise  also  scales  with 
acceleration.  The  result  of  capillary  rise  scaling  is  a  reduction  in  the 
rise  height  and  results  in  desaturation  of  the  pore  fluid. 

Desaturation  curves  for  Tyndall  and  Ottawa  sands  are  shown  in 
Figures  3.9(a)  and  3.9(b).  Capillary  rise  heights  at  1  g  for  the  two 
sands  are  2.04  meters  and  1.S3  meters,  respectively.  X  typical  specimen 
height  for  this  investigation  is  0.28  metars,  and  at  1  g  both  sands  have 
the  potential  to  hold  at  least  70  percent  saturation  to  a  height  of  0.28 
meters.  When  the  specimens  are  accelerated  to  18.9  g's,  the  height  of 
capillary  rise  becomes  0.11  meters  and  0.08  meters  for  Tyndall  and  Ottawa 
sands,  respectively.  At  26.3  g's,  the  height  of  capillary  rise  becomes 
0.08  meters  and  0.06  meters  for  Tyndall  and  Ottawa  sands.  Thus,  the 
saturation  of  the  soil  above  these  capillary  heights  is  at  the  residual 
saturation  which  is  approximately  11  and  5  percent  for  Tyndall  and  Ottawa 
sand  respectively  (Section  III.B.l). 

The  instrument  layer  is  located  at  a  height  of  20  cm  within  the 
specimen.  Thus,  at  an  acceleration  of  16.86  g's,  the  height  of  the 
instruments  is  well  above  the  capillary  rise  height  of  0.11  meters  and 
0.06  meters  for  Tyndall  and  Ottawa  sands.  As  a  check  to  monitor  the 
height  of  capillary  rise  for  Tyndall  sand,  water  content  readings  for 
moist  specimens  were  obtained  at  various  depths  immediately  after  testing. 
Plots  of  saturation  versus  height  within  the  specimens  for  Tyndall  sand 
are  shown  in  Figures  C.l(a)  and  (b).  Curves  plotted  for  Tyndall  sand 
compacted  at  17,  35,  and  S3  percent  saturations  follow  a  typical 
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HEIGHT  iN  SPECIMEN  ycm)  HEIGHT  IN  SPECIMEN  (cm) 


Figur*  C.l  Saturation  veraua  haight  within  soil  spaciotan.  (a)  Spacimana 
acceiwratad  at  18.86  g's.  (b)  Spacimans  aeealacatad  at  26.34 
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dsiaturation  path,  and  capillary  risa  heighta  fall  in  tha  range  from  12 
to  17  cm. 


Due  to  exceusive  water  in  the  70  percent  aaturation  epecimene  of 
Tyndall  sand,  a  special  dry  layer  was  placed  at  the  bottom  of  the  bucket. 
This  allowed  for  a  greater  capacity  of  pore  fluid  to  be  stored  in  the 
bottom  of  the  specimen  away  from  the  instrumentation. 


) 


(The  reverse  of  this  page  is  blank) 
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APPBKDIX  D 


PEAK  STRESS  VERSUS  PEAK  PARTICLE  VELOCITY 


Plots  of  poak  Btrass  as  a  function  of  paak  particla  valocity  for 
Tyndall  beach  sand  are  shown  in  Figura  D.l.  Nota  that  from  Equation 
(2.15) 


P^c 


(D.l) 


and  thus,  the  acoustic  impedance  of  Tyndall  beach  sand  is  equal  to  the 
slope  of  the  regression  through  the  data  in  Figure  D.l.  The  mean 
regression  slope  of  all  the  data  form  Figure  D.l  is  529., 5  kPa/m/s. 
Values  of  acoustic  impedance  determined  using  Equations  (5.3)  and  (5.4), 
for  Tyndall  sand  of  various  saturations,  are  presented  in  Table  5.6a. 
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■  0%  SATURATION  +  17%  SATURATION  *  35%  SATURATION 

a  53%  SATURATION  >:  70%  SATURATION 


Figure  D.l  Plot  of  peak  stress  as  a  function  of  peak  particle  velocity. 


